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Creep  and  Strength  Behavior  of  Frozen  Silt 
in  Uniaxial  Compression 


ZHU  YUANLIN  AND  DAVID  L.  CARBEE 


INTRODUCTION 

The  design  of  stable  structures  in  cold  regions  requires  an  understanding  of  the  creep  and 
strength  behavior  of  frozen  soils.  While  many  researchers  have  studied  this  subject  for  a  long 
time,  detailed  studies  are  still  required  to  gain  a  better  understanding  of  the  subject.  The 
presence  of  ice  and  unfrozen  water  in  frozen  soil  makes  its  stress-strain-strength  behavior  to 
be  strongly  time  and  temperature  dependent  because  the  phase  equilibrium  between  ice  and 
unfrozen  water  is  controlled  by  temperature  and  the  stress  state.  This  temperature  and  stress 
dependency  can  be  evaluated  in  the  laboratory  using  either  constant-stress  or  constant-strain- 
rate  tests.  Mellor  (1980)  pointed  out  that  there  is  a  correspondence  between  these  two  test 
methods  for  ice,  making  them  interchangeable.  Ladanyi  (1981)  reported  that  this  was  also 
true  for  frozen  soil. 

The  purposes  of  this  study  were  to  evaluate  the  influences  of  applied  stress  (or  applied 
strain  rate)  and  temperature  on  the  strength  and  creep  behavior  of  frozen  soil  and  to  system¬ 
atically  investigate  the  correspondence  of  results  obtained  by  constant-strain-rate  and 
constant-stress  tests. 

These  tests  were  conducted  on  a  frozen  silt  at  seven  temperatures  ranging  from  -0.5°  to 
-10°C  and  three  nominal  densities  ranging  from  1.08  to  1.40  g/cm\  The  constant-stress 
creep  tests  were  conducted  at  stress  levels  causing  test  durations  from  a  few  minutes  to  more 
than  two  months.  Constant-head  velocity  tests  were  employed  in  lieu  of  constant-strain-rate 
tests  because  of  difficulties  in  the  testing  procedures.  The  resulting  average  strain  rates 
ranged  from  1.1  x  10~‘  to  6.2xl0‘l  s'1. 

This  report  presents  representative  test  results,  the  analysis  of  the  stress-strain-tempera¬ 
ture  relationships,  and  the  determination  of  the  correspondence  of  the  constant-stress  and 
constant-strain-rate  test  results.  Readers  interested  in  the  raw  data  and  graphs  are  referred  to 
Zhu  and  Carbee  (1983)  available  at  the  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory. 


REVIEW  OF  PREVIOUS  WORK 

Since  Tsytovich  published  the  first  paper  reporting  a  study  of  the  mechanical  properties  of 
frozen  soil  in  1930,  numerous  investigations  in  this  field  have  been  published.  In  1937  and 
1952,  Tsytovich  and  his  co-workers  published  two  monographs  in  which  the  fundamental 
principles  of  frozen  soil  mechanics  were  first  set  forth.  In  1952,  the  Arctic  Construction  and 
Frost  Effects  Laboratory  (ACFEL)  published  a  report  summarizing  test  data  obtained  up  to 
that  time,  including  the  results  of  ACFEL  investigations.  In  later  years,  Tsytovich  (1954, 


1958)  and  Vialov  (1959,  1962,  1963)  published  rather  complete  data  on  the  strength  and  de- 
formational  properties  of  both  remolded  and  undisturbed  frozen  soils.  These  publications 
summarized  and  formulated  qualitative  theories  and  empirical  equations  describing  strength 
and  deformation  of  frozen  soil  as  a  function  of  temperature  and  duration  of  load  applica¬ 
tion.  These  results  are  still  widely  used  in  engineering  practice.  Meanwhile,  Sanger  and  Kap- 
lar  (1963)  published  unconfined  compression  creep  test  data  on  a  variety  of  soils  at  various 
temperatures  from  0°to  -7.8^,  and  presented  empirical  equations  relating  creep  strain  and 
strain  rate  to  applied  stress  and  temperature. 

To  apply  the  rate  process  theory  (RPT)  to  the  creep  process  of  frozen  soil,  Andersland  and 
Akili  (1967)  performed  unconfined  compression  creep  tests  on  a  partially  saturated  frozen 
clay  and  arrived  at  an  activation  energy  of  93.1  kcal/mole  for  a  stress  range  of  42.2-56.3 
kg/cm1  and  temperature  range  of-12°to-18cC.  They  presented  an  empirical  equation  for 
predicting  strain  rate  based  on  the  RPT.  Mitchell  et  al.  (1968)  and  others  have  shown  that  the 
temperature  and  stress  dependency  qualitatively  fit  RPT  predictions  very  well.  Goughnour 
and  Andersland  (1968)  published  uncon  fined  compression  strength  and  creep  data  for  ice 
and  Ottawa  sand-ice  samples  with  various  volumetric  ratios.  An  empirical  equation  relating 
creep  rate  to  stress,  temperature,  strain  and  strain  energy  was  fitted  to  the  ice  sample  data. 
Using  this  equation  for  ice  and  the  stress-strain  curves  for  Ottawa  sand-ice  samples,  they 
plotted  creep  curves  for  sand-ice  specimens  by  means  of  stress  factors.  The  stress  factors 
were  related  to  the  percent  of  sand  by  volume.  Sayles  (1968)  and  Sayles  and  Haines  (1974) 
published  unconfined  compression  creep  data  for  frozen  Ottawa  sand,  Manchester  fine 
sand,  Suffield  clay,  Hanover  silt  and  columnar-grained  ice  for  temperatures  ranging  from 
-0.56°  to  -9.45 °C.  Sayles  found  that  Vialov’s  creep  equation  fit  the  test  data  very  well.  He 
also  developed  a  simplified  method  for  predicting  creep  deformation  that  also  provided  a 
good  fit  to  the  test  data. 

Combining  classical  creep  theories  for  metals  and  existing  creep  theories  for  frozen  soil, 
Ladanyi  (1972)  developed  a  macroanalytical  secondary  creep  model  (an  engineering  theory 
of  creep)  for  frozen  soil,  which  is  simpler  than  Vialov’s  primary  creep  model  for  predicting 
long-term  creep  deformation  and  strength  of  frozen  silt.  Moreover,  Vialov  (1973)  established 
a  micromechanistic  theory  of  creep  deformation  and  failure  of  clayey  soil  based  on  an  in¬ 
vestigation  of  the  kinetics  of  changes  in  the  soil  microstructure.  The  proposed  equations  for 
predicting  creep  deformation  and  long-term  strength  also  fit  their  test  data  very  well. 

In  the  past  decade,  CRREL  has  published  many  reports  on  constant-strain-rate  compres¬ 
sion  tests  on  frozen  soils.  Sayles  (1974)  published  constant-strain-rate  tests  on  frozen  silt  and 
clay  at  an  average  strain  rate  of  0.14-0.15  min'1.  Empirical  power-law  equations  were  pre¬ 
sented  to  evaluate  the  peak  (maximum)  compressive  strength  and  the  50%  peak  strength 
modulus  as  a  function  of  temperature.  Haynes  et  al.  (1975)  investigated  the  strain  rate  effect 
on  the  strength  of  frozen  Fairbanks  silt  at  -9.4<G.  They  concluded  that  the  unconfined  com¬ 
pressive  strength  is  very  sensitive  to  strain  rate  and  increased  10  times  over  a  strain  rate  range 
of  2.9  xlO'4  to  2.9  s'1.  They  reported,  however,  that  the  initial  tangent  modulus  is  not  very 
sensitive  to  strain  rate  and  has  a  magnitude  ranging  from  10*  to  2.2  x  105  kg/cnri.  Later, 
Haynes  and  Karalius  (1977)  studied  the  effect  of  temperature  on  the  strength  of  the  silt  at 
machine  speeds  of  0.0423  and  4.23  cm/s.  They  pointed  out  that  both  the  strength  and  the  ini¬ 
tial  tangent  modulus  are  very  sensitive  to  temperature,  and  they  increased  about  one  order  of 
magnitude  as  the  temperature  was  lowered  from  0°to  -56.7“G.  Empirical  equations  relating 
the  peak  strength  with  unfrozen  water  content  were  also  presented.  Sayles  and  Carbee  (1981) 
investigated  the  effect  of  ice  content  and  dry  unit  weight  on  the  strength  of  frozen  Fairbanks 
silt  at  a  strain  rate  of  5  x  10"’  s'1  and  a  temperature  of  -1 ,67qC.  They  found  that  the  relation¬ 
ship  between  the  compressive  stress  at  the  onset  of  fracturing  and  the  dry  unit  weight  was 
better  defined  than  that  between  peak  strength  and  dry  unit  weight,  and  the  ‘‘ice  matrix 
strength”  was  nearly  proportional  to  the  volumetric  ice  content  of  the  soil. 

2 


More  recently,  a  number  of  researchers  evaluated  the  applicability  of  the  rate  process 
theory  for  frozen  soil.  Analyzing  the  kinetic  nature  of  the  long-term  strength  of  frozen  soils. 
Fish  (1980)  reported  that  the  failure  activation  energy  of  frozen  soils  was  relatively  stable;  it 
varied  only  from  12.9  to  19.7  kcal/moie  in  the  temperature  range  from  -0.55°  to  -20X. 
Based  on  this  analysis,  Fish  derived  a  unified  constitutive  equation  that  can  describe  the  en¬ 
tire  creep  curve  from  the  primary  through  the  tertiary  stage.  This  equation  has  a  form  similar 
to  Assur’s  (1980)  creep  model  for  ice.  Martin  et  al.  (1981)  and  Ting  (1981)  published  rather 
complete  creep  data  for  frozen  Manchester  fine  sand  and  ice  at  a  temperature  range  of  - 1 1  ° 
to  -27X1.  The  free  energy  of  activation  was  reported  by  Martin  et  al.  (1981)  to  have  a  value 
of  1 14  kcal/mole  for  an  ice  saturation  of  40%  and  76  kcal/mole  at  100%.  A  “tertiary”  creep 
model  similar  to  Assur’s  model  was  also  developed  by  Fish  (1980)  and  Ting  (1981). 

Zhu  et  al.  (1982)  reported  an  in-situ  investigation  of  creep  of  massive  ground  ice  with  soil 
at  a  temperature  of  -0.4  <C.  A  simple  power  law  was  suggested  for  describing  the  flow  law  of 
warm  ground  ice.  Wu  et  al.  (1982)  published  a  comprehensive  report  summarizing  the  results 
of  laboratory  circular-footing  creep  tests  for  various  frozen  soils.  They  presented  a  number 
of  empirical  equations  relating  steady-state  creep  rate,  time  to  failure  and  failure  strain  to 
applied  stress,  temperature  and  water  content.  A  f  **  primary  creep  model  fit  decaying  creep 
curves  very  well.  A  simple  power-law  strength-loss  equation  was  also  presented. 

By  analyzing  the  creep  and  strength  test  data  on  frozen  sand  (Martin  et  al.  1981,  Ting 
1981),  Ting  et  al.  (1983)  proposed  and  quantified  various  mechanisms  of  strength  of  frozen 
sand. 


SPECIMEN  PREPARATION 
Material 

The  material  used  in  this  investigation  was  a  remolded  silt  from  the  CRREL  experimental 
permafrost  tunnel  at  Fox,  near  Fairbanks,  Alaska.  Its  gradation  curve  is  shown  in  Figure  1 . 


U.S.  Std.  Sieve  No.  Hydrometer 


Figure  1.  Gradation  curve  of  Fairbanks  silt. 
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Some  of  the  physical  properties  of  the  silt  are  as  follows: 


Plastic  limit  34.2% 

Liquid  limit  38.4% 

Organic  content  5.5% 

Specific  gravity  2.680 

Specific  surface  area  35.0  mVg. 

The  silt  is  classified  as  ML  using  the  Unified  Soil  Classification  System.  The  relationship  be¬ 
tween  unfrozen  water  content  and  temperature  for  this  silt  is  shown  in  Appendix  A. 

Molding 

Distilled  water  was  added 
to  40  lb  of  air-dried  Fair¬ 
banks  silt  to  make  an  initial 
water  content  of  12%  by 
weight.  After  storage  over¬ 
night  to  allow  for  moisture 
equilibration,  the  moist  soil 
was  carefully  compacted  to 
the  desired  density  in  a  gang 
mold.  The  mold,  machined 
from  acrylic  plastic  (Fig.  2), 
could  form  nineteen  7-cm- 
diameter  by  19.3-cm-long 
specimens.  A  specially  de¬ 
signed  compactor  was  used 
to  uniformly  tamp  the  silt- 
water  mixture  in  12  layers, 
each  layer  having  a  thick¬ 
ness  of  1.61  cm.  Samples  Figure  2.  Soil  specimen  gang  mold. 


Figure  2.  Soil  specimen  gang  mold. 
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Figure  3.  Set-up  of  saturation  system. 


Figure  4.  Freezing  cabinet  with  free  water  channel  through  bottom. 


were  compacted  to  nominal  dry  densities  of  1.08,  1.20  and  1.40  g/cm’.  These  dry  densities 
are  referred  to  as  low,  medium  and  high  in  the  remainder  of  this  report. 

After  compaction  the  specimens  were  saturated  with  the  deaerated,  distilled  water  under  a 
vacuum  of  73  mm  Hg.  The  saturation  system  set-up  is  shown  in  Figure  3. 

The  specimen-charged  mold  was  then  placed  into  a  freezing  cabinet  (Fig.  4)  and  quickly 
frozen  from  the  top  down  in  an  open  system  with  an  upper  boundary  temperature  of  lower 
than  -30X.  To  minimize  the  formation  of  ice  lenses  during  freezing,  the  samples  were  cooled 
in  the  cabinet  at  OX  overnight  before  freezing.  Examination  of  the  cut  face  of  specimens  re¬ 
vealed  no  visible  ice  lenses  except  for  those  specimens  with  the  lowest  density,  in  which  some 
ice  lenses  can  be  seen,  as  shown  in  Figure  S. 

After  freezing,  the  specimens  were  ejected  from  the  mold.  One  specimen  from  each  mold 
was  cut  into  three  pieces  to  evaluate  any  density  and  water  content  variations.  The  results  for 
the  representative  samples  from  12  molds  are  shown  in  Table  1 .  For  the  specimens  with  me¬ 
dium  and  low  density,  the  distribution  of  dry  density  and  water  content  along  the  specimen 
length  is  quite  uniform.  However,  for  those  specimens  with  high  density  (molds  12  and  13),  a 
significant  amount  of  water  appears  to  have  migrated  towards  the  top  of  the  specimen  dur¬ 
ing  freezing. 

The  remaining  specimens  were  inspected  for  imperfections,  and  then  the  ends  were  care¬ 
fully  trimmed  flat  and  parallel  on  a  lathe  in  a  coldroom  (Fig.  6).  The  nominal  size  of  the 
specimen  after  trimming  was  70  mm  in  diameter  by  1 52  mm  long.  The  average  squareness  of 
the  specimens  ranged  from  0.1  to  0.2  mm.  The  bulk  density  was  determined,  and  each  speci¬ 
men  was  sealed  with  a  piece  of  T-600  Series  membrane  and  two  steel  end  caps.  Before  test¬ 
ing,  all  specimens  were  tempered  at  the  appropriate  testing  temperature  for  at  least  48  hours. 
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Figure  6.  Trimming  a  specimen  on  a  lathe  in  a  cold  room. 
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Tabic  1.  Dry  density  end  water  content  profile  data  for  12  molds. 


Water  content  (%) _  _ Dry  density  ( g/cm ’) 


Mold 

no. 

Top 

<0-5  cm) 

Middle 
(5-10  cm) 

Bottom 
(10-15  cm) 

Top 

(0-5  cm) 

Middle 
(5-10  cm) 

Bottom 
(10-15  cm) 

i 

42.4 

43.6 

43.7 

1.198 

1.178 

1.181 

2 

44.2 

44.9 

45.2 

1.177 

1.167 

1.162 

3 

41.9 

40.2 

40.6 

1.204 

1.235 

1.227 

4 

43.8 

43.9 

43.0 

1.167 

1.173 

1.182 

5 

40.2 

40.2 

41.3 

1.208 

1.212 

1.200 

6 

40.5 

40.7 

1.228 

1.227 

1.227 

7 

41.4 

42.0 

42.6 

1.218 

1.202 

1.187 

8 

41.5 

40.3 

40.8 

1.216 

1.237 

1.229 

9 

44.7 

44.0 

43.3 

1.163 

1.175 

1.186 

10 

50.7 

51.1 

49.8 

1.082 

1.077 

1.092 

12 

32.1 

29.9 

29.3 

1.388 

1.432 

1.446 

13 

32.7 

28.7 

31.1 

1.368 

1.461 

1.411 

TESTING  PROCEDURE  AND  APPARATUS 

The  uniaxial  constant-stress  creep  tests  were  conducted  using  the  constant-stress  test  ap¬ 
paratus  (Fig.  7)  designed  by  Sayies  (1968).  It  allows  the  applied  load  to  increase  proportion¬ 
ally  to  the  increase  in  the  diameter  as  the  specimen  deforms  so  that  the  true  stress  remained 
constant  during  a  test.  The  test  results  showed  that  the  maximum  deviation  of  the  axial  true 
stress  from  the  initial  applied  stess  is  less  than  1%  of  the  initial  stress  during  tests. 

The  major  portion  of  the  uniaxial  constant-strain-rate  test  program  was  conducted  in  a 
coldroom  on  a  screw-driven  Instron  universal  testing  machine  installed  with  a. temperature- 
controlled  chamber  (Fig.  8).  Some  tests  with  the  highest  machine  speed  (50  cm/min)  were 


Figure  7.  Constant-stress  test  apparatus. 
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a.  Temperature-controlled  testing  chamber. 


b.  Operating  panel  and  recording  system. 

Figure  8.  Instron  universal  testing  machine  model  TT-CM-L  with  an  environmental 
chamber. 


performed  on  a  closed-loop  servo-controlled  electrohydraulic  MTS  testing  machine  (Fig.  9). 
Temperature-controlled  cabinets  were  employed  to  maintain  the  desired  test  temperature. 

The  applied  axial  loads  were  measured  with  various  types  of  load  cells  according  to  their 
appropriate  working  ranges.  The  deformations  were  measured  with  Collins  Corp.  direct-cur¬ 
rent  displacement  transducers,  which  have  a  sensitivity  of  2  x  10'5. 

For  the  creep  tests  and  the  constant-strain-rate  tests  conducted  on  the  Instron  universal 
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figure  9.  MTS  machine  and  Ransco  temperature-controlled  cabinets. 
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Figure  10.  Autodata  Nine  data  logger  with  a  Memodyne  tape  recorder, 
model  3765-8BV. 

machine  in  a  coldroom,  testing  temperatures  were  controlled  by  a  thermistor  positioned  in¬ 
side  the  testing  chamber  around  the  specimens;  the  thermistor  actuated  a  temperature  con¬ 
troller  to  supply  heat  upon  demand.  For  the  constant-strain-rate  tests  performed  on  the 
closed-loop  MTS  machine,  test  temperatures  were  maintained  by  the  Ransco  refrigeration 
unit  and  the  temperature-controlled  cabinets.  During  tests,  the  environmental  temperatures 
around  specimens  were  measured  by  a  thermistor  with  a  sensitivity  of  0.01  £'C.  The  observa¬ 
tions  showed  that  temperatures  were  held  well  within  0.05 'C  of  the  desired  values. 

During  a  test,  the  applied  load,  deformation  and  temperature  were  recorded  with  a  data 
logger,  Autodata  Nine  (Fig.  10),  or  a  strip-chart  recorder,  model  83373-30  (Fig.  1 1),  or  both. 
After  testing,  the  samples  were  photographed,  and  the  bulk  densities  and  water  contents 

were  determined. 


Figure  11.  Strip-chart  recorder,  model  83373-30. 


TEST  RESULTS 


The  detailed  test  results  and  substantiating  graphs  for  each  specimen  are  available  in  a 
CRREL  Internal  Report  entitled  “Constant  Stress  and  Constant  Strain  Rate  Compression 
Test  Data  of  Frozen  Fairbanks  Silt”  (Zhu  and  Carbee  1983).  The  complete  raw  data  can  be 
found  on  magnetic  tapes  at  CRREL.  Only  a  few  representative  graphs  and  summary  tables 
of  the  test  results  are  presented  here.  Tables  2  and  3  summarize  the  creep  test  results  for  satu¬ 
rated  and  partially  saturated  samples,  respectively.  The  constant-strain-rate  test  results  for 
saturated  and  partially  saturated  samples  are  shown  in  Tables  4  and  5.  The  detailed  physical 
data  for  each  specimen  can  be  found  in  Appendix  B.  Some  unspecified  quantities  in  these 
tables  are  defined  as  follows. 


Definition  of  strain  and  stress 

Before  discussing  the  test  results,  it  is  appropriate  to  define  stress  and  strain.  True  (or  nat¬ 
ural)  axial  strain  and  true  axial  stress  are  used  in  this  investigation.  They  are  defined  as  fol¬ 
lows: 


_  .  .  .  final  length 

True  strain  f  =  -In  .-  ,  -r 

mitial  length 


■ ,n  w 


where  ec  is  the  conventional  (or  engineering)  strain,  which  equals  AL/Lo,  where  A L  and  L0 
are  the  axial  deformation  and  the  initial  length  of  specimens,  respectively. 


True  stress  a  = 


axial  load 


instantaneous  cross-sectional  area  A 


If  the  volume  did  not  change  during  testing  and  the  samples  maintained  the  shape  of  right 


Table  2.  Swairy  of  creep  teal  reaelts  for  saturated  temples. 


Specimen 

no. 

Applied 

areas 

o 

(kg/cm‘) 

Minimum 
creep  rale 

(S') 

Time  to 
failure 

(i min J 

Failure 

strain 

_ il— 

6-95 

2.7 

Medinas  density 

9  -  -0.5'C 

3.90x10** 

14,000 

0.113 

6-93* 

3.0 

3.47  x  10-’ 

4,623 

0.152 

9-148 

3.0 

6.11x10” 

329 

0.017 

6-98* 

3.5 

1.46x10-* 

1,504 

0.173 

6-89 

4.0 

4.99x10-* 

111 

0.037 

6-88 

6.0 

3.43x10-’ 

20.3 

0.049 

7-87 

10.0 

3.50x10-* 

3.17 

0.077 

6-86 

12.0 

1.58x10-’ 

0.88 

0.095 

2-18 

4.5 

8  -  -1.0‘C 

2.50x10-* 

27,300 

0.082 

6-94* 

5.0 

3.03x10” 

8,429 

0.181 

9-146 

5.0 

6.36x10” 

1,384 

0.062 

2-20 

5.8 

2.04x10-* 

1,013 

0.161 

6-92 

6.0 

2.31x10” 

637 

0.107 

9-147 

6.0 

2.24x10” 

111 

0.021 

6-91 

7.0 

9.84x10” 

70 

0.051 

4-69 

7.0 

9.50x10” 

66.4 

0.044 

1-13 

9.7 

6.30x10” 

15.7 

0.068 

1-15 

9.7 

6.19x10"’ 

14.5 

0.063 

1-14 

14.2 

4.74x10-* 

3.67 

0.112 

1-12 

18.9 

1.57x10” 

1.08 

0.110 

6-90 

20.0 

1.27x10” 

0.92 

0.083 

2-26 

7.6 

9  •  -2-0°C 

1.80x10” 

29,580 

0.087 

6-97 

7.9 

4.50x10” 

15,470 

0.090 

4-70 

8.5 

8.00x10” 

1,156 

0.069 

2-23 

9.5 

4.28x10-* 

384 

0.112 

6-96 

11.0 

1.04x10” 

115 

0.090 

2-25 

14.4 

6.27x10-’ 

23.8 

0.108 

2-24 

18.9 

2.48x10” 

6.17 

0.109 

3-62 

20.0 

2.97x10” 

4.25 

0.088 

4-75 

24.8 

8.17x10” 

1.67 

0.091 

9-140 

11.0 

9  -  -J.O’C 

3.90x10-* 

5,800 

0.038 

6-100 

11.5 

2.97x10-’ 

1,733 

0.059 

4-80 

12.0 

1.29x10” 

636 

0.075 

4-83 

12.0 

1.13x10” 

573 

0.061 

6-99 

13.0 

2.24x10” 

243 

0.048 

4-77 

15.0 

1.65x10” 

51.5 

0.070 

4-72 

20.0 

1.11x10” 

11.8 

0.094 

4-81 

20.0 

1.18x10” 

11.0 

0.097 

4-71 

25.0 

3.26x10” 

4.5 

0.015 

4-73 

30.0 

8.01  x  10” 

1.83 

0.104 

7-110 

17.0 

•  -  -5.0*C 

1.50x10-’ 

3,864 

0.094 

1-8 

18.0 

8.54x10-’ 

1,064 

0.093 

1-7 

19.1 

1.10x10” 

680 

0.101 

7-102 

20.0 

7.35x10” 

180 

0.092 

*  Failure  did  not  occur  during  test.  Values  shown  for  "Tine  to 
failure”  and  "Failure  strain”  are  for  the  time  and  strain  when 
the  test  was  terminated. 


Specimen 

no. 

Applied 

stress 

a 

( kg/cm V 

Minimum 
creep  rate 

<rt 

Time  to 
failure 
« « 

(min) 

Failure 

strain 

H 

10-155 

20.0 

8.48x10” 

108 

0.070 

7-105 

22.0 

1.81x10*’ 

63.7 

0.090 

1-10 

23.6 

3.67x10-’ 

36.0 

0.107 

1-9 

24.4 

3.82x10-’ 

36.5 

0.102 

1-2 

28.6 

1.32x10” 

11.8 

0.109 

1-4 

28.6 

1.34x10” 

10.7 

0.097 

1-3 

28.6 

1.40x10-* 

11.8 

0.120 

7-103 

30.0 

1.32x10” 

9.83 

0.099 

1-6 

37.7 

6.27x10” 

2.42 

0.098 

1-5 

38.2 

7.34x10” 

2.0 

0.091 

6-101 

40.0 

6.20x10” 

9  -  -7.0*C 

2.0 

0.069 

7-11 

22.0 

3.30x10” 

14,860 

0.074 

7-112 

23.7 

8.62  xlO” 

1.294 

0.097 

7-106 

25.0 

2.11x10” 

454 

0.063 

7-109 

25.0 

2.16x10” 

479 

0.067 

7-107 

30.0 

2.20  xlO” 

58 

0.094 

11-171 1 

35.0 

8.03x10” 

11.2 

0.057 

7-104 

40.0 

2.08x10” 

6.5 

0.100 

ll-I72t 

40.0 

2.60x10” 

9  -  -10.ITC 

3.5 

0.060 

7-114 

32.0 

4.06  xlO” 

1,584 

0.070 

7-113 

35.0 

2.44x10” 

471 

0.103 

7-115 

37.0 

5.04x10” 

224 

0.105 

7-116 

40.0 

l.J9x>l0” 

52.7 

0.071 

il-17St 

50.0 

2.22x10” 

5.20 

0.078 

11-179* 

60.0 

3.60x10” 

1.70 

0.047 

ll-180t 

60.0 

5.14x10"* 

Lew  density 

9  -  -0.5'C 

1.73 

0.054 

10-154 

3.5 

1.48x10” 

9  -  -2.0'C 

80 

10-168 

7.0 

9.00x10” 

974 

0a 

10-164 

7.5 

1.11x10-’ 

453 

0.0 

10-161 

8.0 

2.67x10-’ 

340 

O.OOfce 

10-158 

8.5 

1.10x10” 

168 

0.0160 

10-163 

11.0 

2.98x10” 

80 

0.0206 

10-160 

15.0 

1.80x10-’ 

3.33 

0.0061 

10-157 

16.5 

4.70x10-’ 

2.58 

0.0108 

10-159 

20.0 

1.44x10” 

High  density 

9  -  -2.0  °C 

0.75 

0.0093 

12-202 

10.0 

4.98x10” 

22,940 

0.0863 

12-197 

11.0 

5.00x10-* 

8,478 

0.1030 

12-200 

12.0 

2.19x10” 

460 

0.161 

12-196 

13.0 

7.94x10” 

207 

0.171 

12-195 

15.0 

1 .82  x  10” 

93 

0.168 

12-186 

20.0 

2.69x10** 

6.33 

0.189 

12-199 

25.0 

7.30x10-* 

3.25 

0.191 

t  Tested  on  a  servo-controlled  MTS  machine. 
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Table  3.  Seminary  of  creep  teat  resaHs  for  partially  aatarated  samples  with 
median  density. 


Specimen 

no. 

Saturation 

degree 

*/ 

<%) 

Temperature 

9 

rc> 

Applied 

stress 

a 

(kg/ cm2) 

Minimum 
creep  rate 

ter; 

d  ) 

Time  to 
failure 

u 

(min) 

Failure 

strain 

‘t 

5-SJ 

95.7 

-0.5 

4.0 

6.08x10-* 

103 

0.045 

5-82 

95.3 

-0.5 

6.0 

1.94x10-’ 

22.6 

0.034 

5-84 

95.2 

-0.5 

10.0 

4.74x10-* 

2.17 

0.078 

5-65 

95.5 

-1.0 

8.0 

7.31  x  10"* 

48.3 

0.027 

5-66 

95.3 

-1.0 

9.5 

8.42x10-’ 

15.1 

0.090 

5-64 

95.9 

-1.0 

15.0 

4.69x10- 

3.0 

0.097 

5-63 

95.5 

-1.0 

20.0 

1.61x10’ 

0.83 

0.091 

5-68 

96.2 

-2.0 

8.5 

7.72x10-’ 

598 

0.043 

5-78 

96.3 

-2.0 

10.0 

5.24x10- 

74.0 

0.034 

5-79 

95.7 

-2.0 

15.0 

9.86x10-’ 

8.70 

0.063 

5-67 

96.1 

-2.0 

20.0 

4.10x10-* 

3.17 

0.091 

Table  4.  Sanmary  of  constaat-s  train -rate  compression  test  results  for  saturated  samples. 


Specimen 

no. 

Avenge 
strain  rate • 

e 

(S') 

Peak 

strength 

am 

fkg/cm’J 

Strain 

at  «m 

*/ 

Time  to 
failure 

U 

(min) 

Initial 

yield 

strength 

0- 

(kg/cm2) 

Strain 

at  Oy 

b 

Initial 

tangent 

modulus 

E, 

(kg/cm1) 

50% 

strength 

modulus 

E, 

(kg/cm2) 

14-215* 

6.20x10- 

41.0 

0.0614 

»  »  -a - »^_ 

NMIUI  flCMiy 

9  m  -0J*C 

0.017 

l4-216t 

6.20x10- 

38.2 

0.0612 

0.017 

— 

— 

— 

— 

9-149 

5.85x10- 

21.0 

0.0928 

0.28 

15.4 

0.0058 

3360 

3360 

9-151 

1.06x10- 

14.7 

0.0646 

1.02 

13.1 

0.0059 

2661 

2550 

9-152 

1.14x10- 

13.1 

0.0764 

1.12 

11.1 

0.0058 

— 

2690 

mss 

1.00x10-* 

7.2 

0.0093 

1.55 

6.4 

0.0036 

4017 

2101 

11-184 

9.23  x  10- 

4.1 

0.0144 

26.0 

3.2 

0.0021 

3130 

2005 

14-213* 

6.40  x  10- 

49.8 

0.0774 

9  -  -1.0'C 

0.020 

14-214t 

6.20x10-’ 

58.0 

0.0615 

0.017 

— 

— 

— 

— 

9-150 

5.73  x  10- 

28.6 

0.0757 

0.22 

20.7 

0.0033 

3900 

3900 

2-21 

1.10x10-’ 

21.1 

0.1150 

1.75 

12.8 

0.0046 

4710 

2776 

2-31 

1.12x10- 

18.3 

0.1399 

2.08 

— 

— 

3420 

2790 

2-32 

1.10x10-’ 

18.4 

0.1377 

2.08 

12.5 

0.0044 

4250 

2500 

2-29 

1.12x10- 

9.9 

0.1140 

17.0 

8.3 

0.0037 

7500 

3570 

2-30 

1.11x10- 

10.2 

0.1261 

19.0 

8.4 

0.0041 

4000 

2800 

2-33 

1.00x10-’ 

5.4 

0.0219 

38.0 

4.4 

0.0027 

3978 

1786 

2-28 

1.06x10-’ 

6.1 

0.0323 

51.0 

4.7 

0.0028 

2152 

2152 

3-59 

1.01x10-’ 

5.9 

0.0207 

34.0 

4.6 

0.0026 

3978 

1591 

9-145 

1.06x10- 

4.4 

0.0185 

290 

3.4 

0  0021 

5943 

1453 

14-208 1 

6.20x10-’ 

68.5 

0.0457 

9  -  -2.0*C 

0.013 

14-210* 

6.20x10-’ 

72.7 

0.0460 

0.013 

— 

— 

— 

— 

14-207* 

5.85  x  10-’ 

42.3 

0.0772 

0.22 

35.0 

0.0060 

— 

6330 

14-211* 

5.85x10" 

34.9 

0.0770 

0.22 

27.8 

0.0060 

— 

5100 

3-38 

1.15x10-’ 

26.5 

0.1088 

1.58 

20.6 

0.0056 

5333 

4766 

4-51 

1.07x10-’ 

24.6 

0.0854 

1.33 

18.4 

0.0055 

7130 

5100 

*  Computed  by  i  «  <,/(/„  x  60). 
t  Tested  on  an  MTS  machine. 


Table  4  (coat’d).  Summary  of  co ns tant-s train-rate  compression  test  results  for  saturated  sam 
pies. 


Specimen 

no. 

Average 
strain  rate • 
i 

(s') 

Peak 

strength 

°m 

(kg/ cm') 

Strain 
at  am 

V 

Time  to 
failure 

( min ) 

Initial 

yield 

strength 

°r 

(kg/cm'J 

Strain 

at  Oy 

b 

Initial 

tangent 

modulus 

E, 

(kg/cm') 

50 46 
strength 
modulus 
E, 

(kg/cm') 

3-36 

1.1 1  x  10— 

14.7 

0.1271 

19.0 

10.1 

0.0036 

5258 

2285 

3-40 

1.14x10-* 

15.0 

0.1219 

20.0 

11.3 

0.0039 

7450 

3200 

4-50 

1.13x10-* 

14.9 

0.1208 

20.0 

11.5 

0.0039 

6135 

3200 

2-34 

1.12x10-* 

9.7 

0.0768 

114 

5.7 

0.0026 

3635 

2737 

3-42 

1.11x10- 

10.3 

0.0898 

135 

7.6 

0.0031 

5000 

2455 

3-41 

1.11x10-* 

8.0 

0.0817 

1227 

5.0 

0.0024 

3200 

1650 

3-49 

1.03x10-* 

7.3 

0.0383 

619 

4.3 

0.0023 

2800 

1520 

8  -  -3.0°C 

14-219t 

6.10x10-’ 

93.3 

0.0301 

0.008 

— 

— 

— 

— 

14-220t 

6.10x10-’ 

89.9 

0.0271 

0.008 

— 

— 

— 

— 

9-142 

5.63x10-' 

48.3 

0.0563 

0.17 

39.1 

0.0062 

10560 

7786 

9-143 

1.12x10-’ 

33.7 

0.0752 

1.12 

28.5 

0.0052 

10300 

6073 

3-44 

1.12x10’ 

32.9 

0.1066 

1.58 

23.9 

0.0044 

6490 

4900 

3-43 

1.13x10* 

19.5 

0.1080 

16.0 

14.8 

0.0036 

8750 

4114 

3-47 

1.12x10-* 

20.0 

0.0940 

14.0 

14.0 

0.0037 

6250 

3509 

3-45 

1. 10x10-’ 

13.5 

0.0759 

115 

9.6 

0.0029 

4500 

3395 

3-46 

1.11x10-’ 

13.4 

0.0853 

128 

10.3 

0.0029 

11970 

3248 

3-48 

1.03x10- 

11.0 

0.0239 

387 

9.1 

0.0031 

6875 

3458 

9-144 

1.08x10-* 

11.1 

0.0240 

369 

8.9 

0.0028 

7632 

3500 

8  -  -5.0‘C 

14-217t 

6.00x10-’ 

121.2 

0.0149 

0.004 

— 

_ 

— 

— 

14-218t 

6.00  x  10"’ 

122.5 

0.0149 

0.004 

— 

— 

— 

— 

9-138 

5.62x10-’ 

58.7 

0.0438 

0.13 

— 

— 

— 

— 

8-130 

1.11x10-’ 

37.7 

0.0613 

1.22 

30.5 

0.0041 

11360 

8545 

8-131 

1.12x10- 

41.8 

0.0872 

1.30 

33.3 

0.0041 

12000 

8840 

8-132 

1.15x10* 

25.8 

0.1245 

18.1 

19.5 

0.0031 

14540 

5803 

8-133 

1.06x10-’ 

19.4 

0.0325 

51.2 

15.1 

0.0027 

11400 

5134 

8-134 

1.13x10- 

20.1 

0.0853 

126 

15.5 

0.0029 

10500 

4090 

8-133 

1.07x10- 

17.3 

0.0458 

711 

13.1 

0.0029 

12280 

4518 

8  -  -7.0’C 

14-222t 

6.00x10- 

138.8 

0.0180 

0.005 

_ 

_ 

_ 

_ 

14-223 1 

6.00x10-’ 

135.8 

0.0180 

0.005 

— 

— 

— 

— 

9-139 

5.57x10-’ 

81.2 

0.0334 

0.10 

61.5 

0.0050 

22390 

22390 

8-124 

1.11x10- 

52.0 

0.0651 

1.28 

40.2 

0.0047 

14200 

9430 

8-125 

1.11x10-’ 

52.4 

0.0798 

1.20 

41.9 

0.0044 

16800 

11800 

8-126 

1.05x10- 

33.9 

0.0450 

7.13 

30.1 

0.0039 

15410 

10590 

8-127 

1.15x10- 

23.3 

0.0925 

134 

19.2 

0.0030 

12600 

7500 

8-128 

1.04x10- 

25.2 

0.0446 

71.3 

19.4 

0.0030 

17000 

6070 

8  -  -70.0°C 

14-221 1 

6.00x10-' 

176.4 

0.0210 

0.006 

_ 

_ 

_ 

_ 

U-177t 

5.56x10-’ 

100.3 

0.0222 

0.07 

— 

— 

— 

— 

l!-178t 

5.15x10-’ 

103.9 

0.0206 

0.07 

82.9 

0.0055 

— 

— 

9-137 

5.57x10-’ 

100.1 

0.0334 

0.10 

77.0 

0.0047 

— 

17700 

8-120 

1.11x10-’ 

70.2 

0.0872 

1.31 

54.1 

0.0044 

22600 

14500 

8-121 

1.23x10-’ 

72.5 

0.0971 

1.32 

58.8 

0.0041 

21200 

12260 

7-118 

1.15x10- 

49.1 

0.1208 

17.5 

40.1 

0.0038 

27200 

15200 

8-119 

1.13x10- 

38.8 

0.0993 

147 

30.0 

0.0031 

26910 

12860 

8-122 

1.12x10-’ 

38.6 

0.0872 

130 

22.2 

0.0024 

10000 

9400 

8-123 

1.09x10- 

33.1 

0.0651 

998 

17.8 

0.0024 

9000 

6210 

•  Computed  by  i  =  e,/(/mx60). 
t  Tested  on  an  MTS  machine. 
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Tabic  4  (coat’d).  Summary  of  constant-strain-rate  compression  test  results  for  saturated  sam- 


pies. 

Average 
strain  rale' 
Specimen  i 

no.  (s') 

Peak 

•  strength 
a  m 

(kg/cm1) 

Strain 

o< 

Time  to 
failure 
t » 

(min) 

Initial 

yield 

strength 

By 

(kg/cm') 

Strain 

at  By 

_ It _ 

Initial 

tangent 

modulus 

E  i 

(kg/cm') 

SO % 
strength 
modulus 

E, 

( kg/cm V 

11-185 

5.85x10-* 

41.5 

0.0246 

Low  density 

0  =  -2.0°C 

0.07 

36.7 

0.0057 

6600 

10-165 

5.58x10-* 

37.5 

0.0106 

0.04 

32.0 

0.0050 

_ 

7100 

10-166 

1.00x10-' 

29.2 

0.0093 

0.16 

23.9 

0.0041 

_ 

7600 

10-167 

1.00x10- 

15.8 

0.0082 

1.50 

13.6 

0.0026 

11700 

5030 

10-169 

1.01  xlO- 

9.3 

0.0146 

24.0 

7.0 

0.0018 

21200 

3214 

10-162 

8.1x10-’ 

8.0 

0.0102 

210 

6.2 

0.0019 

5350 

3125 

14-203 1 

6.22x10* 

56.0 

0.0933 

High  density 

0  =  -2.0°C 

0.025 

14-204 1 

6.09x10-* 

50.7 

0.1096 

0.029 

_ . 

_ 

_ 

_ 

14-205 1 

6.55  x  10-* 

37.7 

0.1965 

0.50 

21.6 

0.0060 

_ 

_ 

12-188 

6.14x10'* 

39.3 

0.1842 

0.50 

20.7 

0.0060 

6670 

953 

12-189 

1.24x10-’ 

25.5 

0.2046 

2.75 

10.5 

0.0040 

_ 

_ 

12-190 

1.24x10- 

27.8 

0.2146 

2.88 

13.4 

0.0049 

3550 

980 

12-191 

1.24x10-* 

17.6 

0.2243 

30.0 

8.0 

0.0036 

5000 

732 

12-192 

1 .23  x  10"' 

13.1 

0.1863 

253 

5.4 

0.0028 

6000 

475 

12-193 

1.22x10-’ 

12.3 

0.1692 

231 

5.0 

0.0028 

4750 

300 

12-194 

1.19x10" 

11.6 

C  1504 

2098 

4.7 

0.0032 

4006 

396 

12-198 

1.18x10- 

11.8 

0.1507 

2126 

5.6 

0.0037 

5000 

438 

•  Computed  by  i  »  <t/(/mx60). 
t  Tested  on  an  MTS  machine. 

Table  S.  Summary  of  constant-strain-rate  compression  test  results  for  partially  saturated  samples  with  medium 

density. 

Saturation 

degree  Temperature 

Average 

strain 

rate 

Peak 

strength  Strain 

Time  to 
failure 

Initial 

yield 

strength  Strain 

Initial 

tangent. 

modulus 

50% 

strength 

modulus 

Specimen 

»/ 

0 

( 

«m  at 

By 

at  By 

E, 

E, 

no. 

(%) 

rcj 

(sr') 

( kg/cm ')  ir 

(min) 

(kg/cm')  t. 

(kg/cm') 

(kg/cm1) 

5-57 

96.0 

-1 

5-58 

95.7 

-1 

5-54 

96.0 

-4 

5-56 

96.2 

”4 

5-52 

96.2 

5-53 

98.1 

circular  cylinders,  the  true  axial  stress  can  be  calculated  by 


a  =  P(1-€c)//4, 


where  A,  is  the  initial  cross-sectional  area  of  samples. 


Definition  of  creep  failure 

Figure  12,  a  typical  log/  vs  log/  curve  for  a  relatively  high  stress,  clearly  shows  that  the  ax¬ 
ial  strain  rate  t  changes  with  time  /  and  attains  a  minimum  creep  rate  at  a  definite  time.  Fol- 


7700 

4150 

4930 

3600 

?* 

& 


j* 

s!. 


Figure  12.  Typical  curve  of  logl  vs  logl  under  relatively  high 
stress. 


Figure  13.  Typical  curve  of  logl  vs  logl  under  relative¬ 
ly  low  stress  for  medium-density  samples. 


lowing  the  definition  of  creep  failure  by  many  researchers  (Fish  1979,  Assur  1980,  Martin  et 
al.  1981,  Mellor  and  Cole  1982),  we  define  creep  failure  as  the  point  at  which  the  creep  rate 
reaches  its  minimum  value  #m.  Correspondingly,  the  time  to  the  minimum  is  the  time  to  creep 
failure  tm,  and  the  strain  at  the  minimum  is  the  failure  strain  tf.  These  three  quantities— lm, 
tm  and  (f — are  designated  as  the  creep  failure  parameters. 

Because  the  strain  rate  curves  for  the  medium-density  samples  tested  at  relatively  low 
stress  (Fig.  13)  fluctuated  frequently  within  the  so-called  secondary  creep  stage,  it  was  diffi¬ 
cult  to  identify  a  minimum  strain  rate  on  this  type  of  curve.  To  determine  the  minimum 
strain  rate  on  this  type  of  curve  in  a  consistent  manner,  we  smoothed  out  the  “spikes”  in  the 
log#  vs  log#  curves  by  averaging  strain  rates  over  larger  time  intervals,  as  shown  in  Figure  13. 

Definition  of  failure  in  constant-strain-rate  tests 

Figure  14  shows  a  typical  set  of  stress-strain  curves  for  various  strain  rates  at  a  tempera¬ 
ture  of  -3^.  Each  curve  has  a  definite  maximum.  Failure  in  uniaxial  compression  is  thus 
defined  as  the  point  at  which  the  true  axial  stress  attains  its  maximum  value  om.  Similarly, 
the  time  to  the  peak  is  called  the  time  to  failure  tm,  and  the  strain  at  the  peak  is  the  failure 
strain  «f. 
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Figure  14.  Typical  stress-strain  curves  for  Figure  15.  Stress-strain  curves  shown  in  Figure 
specimens  with  medium  density  under  various  t  14  expanded  for  the  the  initial  6V»  strain, 
at  -3°C. 


Figure  16.  Determination  of  50V*  peak 
strength  modulus. 


1 


Definition  of  initial  yield  strength 

Following  Sayles  and  Carbee  (1981),  we 
define  the  initial  yield  strength  oy,  as  shown 
in  Figure  IS,  as  the  stress  at  which  the  slope 
of  a  stress-strain  curve  starts  to  decrease  per¬ 
ceptibly  from  the  initial  tangent  to  the  curve, 
which  is  at  the  strain  where  the  pore  ice  would 
be  expected  to  fracture.  The  strain  at  this 
stress  is  defined  as  the  initial  yield  strain  ty. 
Although  the  determination  of  ay  is  not  as 
precise  as  one  would  like,  meaningful  data 
can  still  be  obtained  as  long  as  a  consistent 
procedure  is  used  in  determining  the  yield 
strength  for  all  expanded  stress-strain  curves. 

Determination  of  initial  tangent  modulus 
and  50%  peak  strength  modulus 

It  is  well  known  that  the  determination  of 
the  initial  tangent  modulus  of  a  nonelastic 
material  based  on  a  stress-strain  curve  is  dif¬ 
ficult,  particularly  when  the  curve  is  rounded 
as  shown  in  Figure  16.  As  far  as  we  know, 


,  Strain  Rata  (s'1) 


there  is  no  standard  procedure  for  determining  the  initial  tangent  modulus.  The  reported 
values  of  initial  tangent  modulus  E\  in  Tables  4  and  5  are  determined  graphically  from 
stress-strain  curves.  The  50%  peak  strength  modulus  E\  is  defined  as  the  tangent  modulus  of 
a  point  on  a  stress-strain  curve  at  which  a  —  am/2,  as  is  illustrated  in  Figure  16.  These  values 
are  also  determined  graphically. 


CREEP  BEHAVIOR 

General  nature  of  the  creep  process  and  the  failure  mode 

The  variation  of  strain  rate  in  the  creep  process  can  be  illustrated  with  a  log*  vs  log/  plot. 
Typical  log*  vs  log/ curves  for  various  test  conditions  are  presented  in  Figure  17.  All  of  these 
curves  reveal  the  same  general  law  of  the  creep  process:  When  a  constant  stress  is  applied  to  a 
specimen,  the  creep  rate  continuously  decreases  with  time  (strain  hardening  dominates), 
reaches  a  minimum  value  (point  F  in  these  figures),  and  then  starts  to  accelerate  until  frac¬ 
ture  or  plastic  failure  occurs  (strain  softening  dominates).  That  is,  there  must  be  an  equilibri¬ 
um  (or  critical)  stress-strain  state  between  these  two  opposite  processes— strain  hardening 
and  strain  softening.  At  this  moment,  samples  must  be  in  equilibrium  with  the  external 
stress.  In  this  sense,  it  is  reasonable  to  define  the  minimum  point  on  log*  vs  log/  curves  as  the 
creep  failure.  The  creep  rate  seldom  remains  constant  within  the  secondary  creep  stage. 
Therefore,  the  so-called  steady-state  creep  may  exist  only  at  one  point,  at  least  for  the  inves¬ 
tigating  conditions  reported  here. 

Examining  all  of  the  log*  vs  log/  curves  for  low  stresses,  we  found  that  there  was  a  large 
decrease  in  creep  rate  at  the  onset  of  secondary  creep  for  samples  with  medium  and  high  den¬ 
sities  but  not  for  the  low-density  samples.  This  suggests  that  different  deformation  mechan¬ 
isms  dominate  the  creep  process  of  ice-rich  low-density  samples  than  for  dense  frozen  silt 
samples. 


Figure  17.  Log  l  vs  logt  curves  for  different  conditions. 
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Figure  18.  Typical  Fairbanks  silt  specimens  after  short-term  creep  tests. 


Another  observation  was  the  occurrence  of  spikes  in  the  log*  vs  log*  curves  within  the  sec¬ 
ondary  creep  stage  for  medium-density  samples  tested  at  low  stresses  (Fig.  17a-g).  The  spikes 
occurred  over  a  time  of  a  few  days  to  more  than  one  month  even  though  the  stresses  and  tem¬ 
peratures  were  constant.  The  spikes  did  not  occur  for  both  high-  and  low-density  tests  for  the 
same  orders  of  strain  rate  (Fig.  17h,  i).  This  evidence  indicates  that  the  spikes  are  not  arti¬ 
facts  of  the  test  apparatus  but  are  sample  responses.  The  same  phenomena  (both  the  drop  of 
the  creep  rate  and  the  occurrence  of  spikes  in  secondary  creep)  have  been  observed  on  a  fro¬ 
zen  sand  (Martin  et  al.  1981). 

The  constant-stress  creep  tests  showed  that  all  specimens  failed  plastically  under  the  ranges 
of  stress  and  temperature  employed.  No  brittle  failures  were  observed.  However,  in  examin¬ 
ing  the  failed  specimens,  we  found  that  different  failure  features  could  be  distinguished,  even 
though  they  all  appeared  to  fail  plastically.  For  example,  for  the  short-term  creep  test  (less 
than  1  day  to  failure)  all  test  specimens  failed  plastically  without  visible  cracks,  even  for 
strains  greater  than  25%  and  low  temperatures  (Fig.  18).  However,  for  the  long-term  creep 
tests  (greater  than  about  1  day  to  failure),  after  large  plastic  strains,  x -shaped  shearing 
cracks  were  visible  on  the  sample  surface  as  shown  in  Figure  19,  especially  for  dense  frozen 
silt.  The  same  phenomenon  was  observed  by  Vialov  (1962)  for  a  frozen  dense  clayey  soil.  We 
call  this  plastic  failure  with  cracks. 

Minimum  creep  rate 

In  the  following  section,  we  will  concentrate  on  the  quantitative  evaluation  of  the  mini¬ 
mum  creep  rate  *m  and  time  to  creep  failure  tm  as  a  function  of  stress  and  temperature,  since 
they  are  the  basic  parameters  chosen  in  this  study  for  predicting  creep  deformation  and  long¬ 
term  strength  of  frozen  soils.  We  found  from  this  study  that  the  minimum  creep  rate  strongly 
depends  upon  applied  stress,  temperature  and  dry  density  of  the  soil. 
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Figure  19.  Typical  Fairbanks  silt  specimens  after  long-term  creep  tests,  showing  x  -shaped  shear 
cracks. 
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Stress  dependence 

The  minimum  creep  rate  *m  as  a  function 
of  stress  a  at  various  test  temperatures  for 
medium  dry  density  is  plotted  in  Figure  20 
on  a  log-log  scale.  The  curves  in  this  figure 
are  not  straight  lines,  but  they  can  be  consid¬ 
ered  a  set  of  quasi-parallel  curves.  The  slopes 
of  these  curves  perceptibly  increase  at  a  cer¬ 
tain  minimum  strain  rate  as  stresses  decrease. 
This  can  be  seen  more  clearly  in  plots  of 
logim  vs  1/a  (Fig.  21).  This  figure  shows  a 
family  of  bilinear  curves  that  all  deflect  at 
about  the  same  minimum  strain  rate  (10~‘ 
s'1).  For  convenience  in  describing  these 
data,  we  define  a  critical  creep  rate  ic  that 
can  be  considered  the  minimum  strain  rate 
for  distinguishing  between  two  types  of  creep 
of  frozen  soil:  short-term  creep,  which  has  a 
minimum  creep  rate  greater  than  tc,  and 
long-term  creep,  which  has  a  minimum  creep 
rate  less  than  ic.  The  stresses  corresponding 
to  tc  for  various  temperatures  are  defined  as 
the  critical  creep  stresses.  Moreover,  these 
two  types  of  creep  are  governed  by  different 
deformation  mechanisms,  which  will  be  discussed  later.  In  addition,  Figure  21  shows  that  at 
high  strain  rates  all  of  the  curves  converge  at  a  common  point.  Thus,  the  constitutive  equa¬ 
tions  for  the  frozen  silt  with  medium  density  can  be  described  by  the  following  exponential 
equations. 


Figure  20.  Log-log  plot  of  tm  vs  a  for  vari¬ 
ous  temperatures. 


-jr,  Reciprocal  of  Stress  (kg/cm*) 

Figure  21.  Log  lm  vs  l/o  curves  for  various  temperatures. 


For  short-term  creep  ( a  >  ac), 

*m  =  *.exp|^  -  (4) 

where  /.  and  a,  are  reference  values  independent  of 
temperature  and  are  used  for  convenience  in  present¬ 
ing  this  data.  By  linear  regression  analysis,  we  obtained 
t,  =  8.84  x  10~5  s'1  and  a.  =  71.4  kg/cm1  for  the  tested 
silt  with  dry  density  rd  ranging  from  1 . 1 8  to  1 .23  g/cm5.  The  parameter  k  is  plotted  as  a  func¬ 
tion  of  temperature  in  Figure  22  and  can  be  determined  by 

k  =  53. 1(0/0.)° 72  for  -0.5°  a  0  a  -2°C 
and  (5) 

k  =  42.4(0/0.)’ 02  for  -2°  a  0  z  -IX: 

where  0  is  the  test  temperature  in  °C,  and  0o  is  a  reference  temperature  taken  as  -1  *C. 

For  long-term  creep  (a  s  aj, 


where  zc  is  the  critical  creep  rate,  equal  to  10'*  s'1,  and  ac  is  the  critical  creep  strength  defined 
above.  The  variation  of  ac  with  temperature  is  illustrated  in  Table  6  and  Figure  23  and  can  be 
expressed  as 

ac  =  5.2(0/0,)07*  (7) 

where  oc  is  in  kg/cm1. 

The  parameter  k '  in  eq  6  as  a  function  of  temperature  is  also  plotted  in  Figure  22  and  can 
be  determined  by 


Table  6.  Values  of  ot  iu  eq  6. 


9 

ra 

Oc 

(kg/cm') 

9 

ra 

Oc 

( kg/cm V 

-0.5 

3.2 

-5 

18.4 

-1 

5.2 

-7 

23.5 

-2 

8.5 

-10 

32.5 

-3 

11.7 

k'  =  134(0/0,)"*. 


23 


(8) 


Figure  24.  Power  n  as  a  function  of  temper¬ 
ature. 


Figure  25.  Plot  of  log!  vs  logo  for  Callov- 
ian  silty  sandy  loam.  (After  Vialov  1962, 
Table  7.) 


The log/m  vs  logo  curves  in  Figure  20  can  also  be  represented  approximately  as  a  set  of  bi¬ 
linear  curves,  which  deflect  at  the  strain  rate  tc  of  10'*  s'1.  Thus,  assuming  straight  lines,  the 
stress  dependence  of  /m  for  medium  density  can  be  also  expressed  by  a  power  equation: 

tm  =  tc(a/OQ)n  (9) 

where  the  exponent  n  depends  upon  temperature.  Its  variation  with  temperature  for  both 
long-term  and  short-term  creep  is  shown  in  Figure  24  and  can  be  evaluated  by 

n  m  5.59(0/0,)°  223  (10) 

for  short-term  creep  and 

n  =  27.36(0/0,)°  233  (11) 

for  long-term  creep. 

A  similar  power  function  was  reported  by  Ladanyi  (1972)  to  evaluate  secondary  creep 
rates  as  a  function  of  stress,  based  on  work  in  metals  by  Hult  (1966).  However,  the  param¬ 
eter  oc  referred  to  a  proof  stress  in  his  equation  is  not  defined  as  discussed  here. 

Vialov  (1962)  also  observed  a  consistent  deflection  on  log/  vs  logo  curves  at  a  strain  rate  of 
near  10"‘  s'1  for  a  comparable  soil  (Fig.  25).  In  addition,  the  slope  of  the  log/  vs  loga  curves 
for  short-term  creep  from  Vialov’s  data  is  also  very  close  to  that  from  our  results.  However, 
there  is  a  substantial  difference  between  the  results  for  long-term  creep;  the  slope  is  much 
greater  for  our  data  than  for  Vialov's. 

Temperature  dependence 

In  recent  years,  attempts  have  been  made  to  apply  the  rate  process  theory  (RPT)  to  the 
creep  of  frozen  soil  (Akili  1966,  1970,  Andersland  and  Akili  1967,  Goughnour  1967,  Anders- 
land  and  Douglas  1970,  Vialov  1973,  Mitchell  1976,  Fish  1980,  Martin  et  al.  1981,  Ting 
1981).  The  different  individuals  have  presented  various  formulas  relating  creep  rate  to  tem- 
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Vt  (s'1  K"1) 


as  a  function  of  stress. 


perature  and  stress  based  on  RPT.  To  evaluate  the  applicability  of  this  theory  to  our  test  re¬ 
sults,  the  short-term  creep  test  data  for  medium  density  were  plotted  in  the  log  (tm/T)  vs  1  IT 
plane  as  shown  in  Figure  26.  As  expected,  all  of  the  isostress  curves  are  straight  lines  within  a 
temperature  range  of  -0.5°  to  -7°C.  All  these  isostress  curves  appear  to  converge  at  a  com¬ 
mon  point.  Thus,  we  can  present  a  modified  RPT  expression  for  frozen  silt  as  follows: 


or 


,  XKT 
tm  =  —  exp 


lm  =  AXT  exp 


[  Utt 

’l 

1\1 

T 

T.)\ 

,1  u'i 

i  1 

1|] 

1  R 1 

\T 

T.)  | 

02) 

(13) 


where  X 
K 
T 
h 

<4 

R 

T, 

A, 


constant  (for  frozen  silt  with  medium  density,  X  =  2.4  x  10*'5) 

Boltzman’s  constant  (1.38x10*”  J/K) 
absolute  temperature  (K) 

Planck’s  constant  (6.624  x  10"’*  Js) 

activation  energy  (kcal/mole) 

gas  constant  (1.987  x  10*J  kcal/mole  K) 

characteristic  temperature  (for  frozen  silt  with  medium  density,  T,  =  273.15  K) 
XK/h  =  constant  (for  frozen  silt  with  medium  density,  A,  -  5.0  x  10*’  s'1  K*’). 


To  calculate  Ue,  rearranging  eq  13,  we  have 


where  A  =  -logA,  =  4.3. 

The  l/e  values  calculated  from  the  creep  test  data  as  a  function  of  stress  are  shown  in  Fig¬ 
ure  27;  there  is  a  close  relationship  between  activation  energy  and  stress.  The  activation 
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energy,  then,  is  a  useful  characteristic  by 
which  creep  test  data  can  be  normalized. 
Explicitly,  Ue  strongly  depends  on  stress, 
especially  at  the  lower  stress  range  tested 
in  this  study.  Therefore,  the  values  of  Ue 
without  specifying  stress  have  no  signifi¬ 
cance.  The  explicit  nonlinearity  between 
Ut  and  a  indicates  that  the  relation  Ut  = 
A F-  &a,  which  was  derived  based  on  an 
approximation  of  RPT  (i.e.  the  energy 
supplied  by  the  action  of  external  stress  is 
greater  than  the  thermal  energy)  and  ex¬ 
perimentally  verified  for  a  limited  range 
of  values  by  a  number  of  investigators 
(e.g.  Mitchell  et  al.  1968,  Ting  1981),  is 
not  valid  over  a  wide  range  of  stresses  dis¬ 
cussed.  The  reason  for  the  discrepancy 
may  be  that  this  approximation  is  not  true 
for  low  stress  or  that  more  than  one  mech¬ 
anism  could  be  operating  simultaneously 
in  the  creep  process  at  low  stress  levels. 
By  replotting  Figure  27  for  the  coordinates  of  log  Ut  vs  log(c/cr,  + 1)  as  shown  in  Figure  28, 
we  found  that  experimental  energy  Ue  as  a  function  of  stress  a  can  be  well  described  by 

Ut  -  AF(a/a,  +  l)"4  (15> 

where  a,  =  reference  stress,  taken  as  1  kg/cm1 
6  =  constant,  equal  to  1 .25 

AF  =  free  energy,  defined  in  this  study  as  the  apparent  minimum  energy  required  to 
overcome  the  energy  barrier  when  the  external  stress  is  zero. 

For  the  frozen  silt  tested,  AFwas  found  to  be  as  high  as  8600  kcal/mole,  which  is  almost  two 
orders  of  magnitude  greater  than  the  free  energy  of  activation  reported  by  some  investigators 
(Andersland  and  Douglas  1970,  Martin  et  al.  1981). 

Following  a  modified  RPT  equation  for  ice,  which  includes  a  power  function  to  describe 
the  stress  dependence  (Glen  1955,  Gold  1973,  Langdon  1973,  Weertman  1973,  Homer  and 
Glen  1978),  we  wrote  another  modified  RPT  equation: 


Figure  28.  Plot  of  log  U,  vs  log(o/o ,  + 1). 
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(16) 


where  X '  -  dimensionless  constant 
a.  =  reference  stress 
n  =  exponent  of  <j/ctc,  defined  in  eq  9 
Ut  =  apparent  activation  energy. 


To  determine  t/a  and  X',  eq  16  is  rearranged  as 
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F/gure  29.  Plot  of  logflm/T(o/oJ°J  vs  1/ T. 


Figure  30.  Plot  of  logac  vs  temperature. 


The  creep  test  data  for  a  z.  ac  was  plotted  in  Figure  29  on  the  basis  of  this  equation.  This  plot 
does  not  follow  a  straight  line  as  expected,  indicating  that  there  are  no  identical  values  of  ap¬ 
parent  energy  l/a  and  the  parameter  A"'  for  the  entire  range  of  testing  temperatures.  In  other 
words,  both  £/a  and  X '  can  be  approximately  considered  as  constants  only  within  a  certain 
range  of  temperature.  For  example,  if  the  curve  in  Figure  29  approximates  a  straight  line 
within  the  temperature  range  of  -3°  to  -10°C,  then  l/a  =  81  kcal/mole  and  X'  =  3.0 x  10“ 
However,  for  higher  temperatures  the  curve  becomes  so  steep  that  it  is  difficult  to  calculate 
the  values  of  t/a  and  X' .  This  means  that  eq  16  does  not  work  for  higher  temperatures. 

Combining  his  creep  law  and  RPT,  Ladanyi  (1972)  developed  a  simple  way  to  determine 
an  apparent  activation  energy: 


(/,  =  2.303x273 *nR 


A  logac 


where  aQ  is  the  stress  corresponding  to  a  given  im,  and  0  is  temperature  (°C).  To  evaluate  l/a 
with  this  equation,  a  plot  of  log<rc  vs  6  was  constructed  from  our  test  data  (Fig.  30).  Clearly 
this  plot  is  nonlinear,  illustrating  that  eq  18  is  appropriate  for  a  certain  range  of  lower  tem¬ 
peratures  but  not  for  high  temperatures.  The  parameter  n  in  eq  18  is  also  temperature  depen¬ 
dent.  Taking  an  average  value  of  n  =  8,  we  computed  Ua  =  386  kcal/mole  for  the  tested 
frozen  silt  within  a  temperature  range  of  -5°  to  -  10°C.  Unfortunately,  we  could  not  compare 
this  computed  value  with  that  shown  by  Ladanyi  (1972),  because  no  temperature  data  are 
available  in  his  report. 


Effect  of  dry  unit  weight 

The  minimum  creep  rate  as  a  function  of  dry  density  for  various  stress  levels  at  -2°C  is  il¬ 
lustrated  in  Figure  31.  Obviously,  dry  density  has  no  significant  effect  on  4m  for  relatively 
high  stresses,  but  it  significantly  affects  for  lower  stresses.  For  example,  im  decreases  by 
three  orders  of  magnitude  as  increases  from  1.07  to  1.40  g/cm’  at  a  stress  level  of  10 
kg/cm1,  while  it  remains  almost  the  same  for  the  same  range  of  at  a  stress  of  20  kg/cm/ 
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Figure  31.  Minimum  creep  rates  as  a  func - 
tion  of  dry  density  for  various  stresses  at 
-2°C. 


Figure  32.  Loglm  vs  1/a  curves  for  the 
specimens  with  high  density  at  -2°C. 


Figure  33.  Logt„  vs  log(o/o,)  curves  for  the 
specimens  with  low  density  at  -2  °C. 


The  effect  of  dry  density  on  creep  behav¬ 
ior  can  also  be  seen  by  comparing  Figures  32 
and  33  with  Figure  21.  Figure  32  is  a  plot  of 
log£m  vs  1/a  for  high-density  samples  at 
-2 °C,  showing  a  curve  similar  to  that  for  me¬ 
dium  density.  The  critical  creep  rate  for  sam¬ 
ples  with  high  density  is  also  about  10~‘  s~', 
and  the  critical  creep  strength  ac  is  about 
11.5  kg/cm 2  at  -2°C.  However,  samples  with 
low  density  possess  a  different  creep  behav¬ 
ior.  As  shown  in  Figure  33,  the  log£m  vs 
Iog(a/at)  curve  for  low  density  is  nearly  a 
straight  line  over  a  wide  range  of  tm,  which  is 
similar  to  the  creep  response  of  ice.  Conse¬ 
quently,  the  creep  law  of  ice-rich  frozen  silt 
can  be  described  by  a  simple  power-law 
equation: 

*m  =  ci(a/ai)n  (19) 

where  Ci  is  the  mimimum  creep  rate  when  a 
=  a,.  For  the  frozen  silt  tested,  c,  = 


1.68  x  10'IS  s'1  and  n  =  6.90  for  0  =  -2°C. 
The  change  in  creep  behavior  due  to  the 
variation  of  density  is  attributed  to  the 
change  in  deformation  mechanism.  According  to  Sayles  (1973)  and  Sayles  and  Carbee 
(1981),  the  stress  resistance  of  a  saturated  frozen  silt  is  composed  of  the  bonding  force  and 
the  frictional  resistance  between  soil  particles  combined  with  that  of  the  fractured  ice  crys¬ 
tals.  For  samples  with  different  densities,  different  components  of  the  stress  resistance  domi- 


nate  the  creep  process.  At  low  density,  most  of  soil  particles  are  suspended  in  an  ice  matrix, 
so  the  deformation  behavior  of  the  ice  matrix  dominates  the  creep  process.  This  explains  why 
the  ice-rich  frozen  silt  has  a  creep  behavior  similar  to  that  of  ice.  F  r  dense  frozen  soil, 
however,  with  the  increase  in  the  number  of  contacts  between  mineral  particles,  the  frictional 
force  between  particles  (Sayles  1973)  becomes  more  important  for  governing  the  creep  proc¬ 
ess.  There  may  exist  a  threshold  for  breaking  the  bonding  force  between  soil  particles.  When 
the  external  shear  stress  acting  on  the  contact  areas  is  high  enough  to  overcome  this  thresh¬ 
old,  the  glide  between  soil  particles  will  take  place.  In  this  case,  the  frictional  resistance  be¬ 
tween  soil  particles,  as  well  as  between  fractured  ice  crystals,  will  dominate  the  creep  rate, 
and  the  resistance  is  not  so  sensitive  to  the  sample  density  for  frozen  silt,  as  is  indicated  by 
the  curves  in  Figure  31  for  higher  stresses.  When  gliding  starts,  plastic  deformation  will  con¬ 
tinuously  grow  until  creep  rupture  occurs,  resulting  in  a  type  of  typical  ductile  failure.  This 
type  of  creep  is  called  glide  creep  in  this  report. 

We  hypothesize  that  when  external  shear  stress  is  less  than  the  threshold,  soil  particles  can¬ 
not  glide  by  each  other.  In  this  case  the  creep  process  is  controlled  by  the  dislocation  of  pore 
ice,  and  the  dislocation  velocity  determines  the  creep  rate.  This  type  of  creep  is  called  disloca¬ 
tion  creep  in  this  report.  For  the  dense  frozen  soil,  the  closely  packed  soil  particles  greatly 
impede  the  movement  of  dislocation  of  the  pore  ice,  so  that  the  secondary  creep  rates  at 
lower  stresses  are  much  smaller,  as  observed.  On  the  basis  of  dislocation  theory,  the  disloca¬ 
tion  is  rapidly  multiplied  along  slip  planes  with  increasing  strain.  The  increase  in  the  number 
of  dislocations  greatly  weakens  the  shear  resistance  of  frozen  soil  on  the  slip  planes,  so  that 
dislocation  creep  usually  leads  to  a  type  of  plastic  failure  with  visible  x  -shaped  cracks,  as 
shown  in  Figure  19,  especially  for  high  densities. 

Based  on  this  interpretation,  the  short-term  creep  is  primarily  controlled  by  glide  creep, 
while  the  long-term  creep  is  governed  by  dislocation  creep.  Also,  this  explanation  gives  the 
critical  creep  strength  ac  a  clear  physical  meaning:  It  is  the  macroanalytical  measure  of  the 
threshold  for  breaking  bonding  force  between  soil  particles.  Theoretically,  the  bonding  force 
increases  with  the  decrease  in  unfrozen  water  content  and  the  distance  between  mineral  parti¬ 
cles.  According  to  this,  the  critical  creep  strength  should  increase  with  decreasing  tempera¬ 
ture  and  increasing  dry  density.  In  the  preceding  discussion,  we  have  shown  this  to  be  true. 

Effect  of  saturation 

The  partially  saturated  samples  discussed  here  refer  to  the  samples  saturated  under  normal 
atmospheric  pressure  at  a  nominal  dry  density  of  1 .2  g/cm\  The  original  purpose  for  making 
these  samples  was  to  examine  the  degree  of  saturation  of  samples  saturated  without  benefit 
of  a  vacuum.  We  expected  the  degree  of  saturation  to  be  considerably  lower  than  100%. 
However,  it  was  higher  than  expected.  As  shown  in  Tables  3  and  5,  the  degree  of  saturation  S 
is  as  high  as  95-%%.  Comparing  the  creep  test  results  of  the  saturated  and  the  partially  sat¬ 
urated  samples  for  the  same  test  conditions  (Tables  2  and  3),  we  found  that  the  minimum 
creep  rates  of  partially  saturated  samples  are  consistently  higher  than  those  of  partially  sam¬ 
ples,  even  though  the  saturation  degree  for  the  former  is  only  3-4%  less  than  that  for  the  lat¬ 
ter.  This  may  be  because  there  is  less  ice  and  therefore  greater  stress  in  the  ice  in  the  partially 
saturated  samples,  so  the  creep  rate  becomes  higher. 

Time  to  creep  failure 

The  test  results  here  and  elsewhere  indicate  that  the  time  to  creep  failure  fm  strongly  de¬ 
pends  on  stress,  temperature  and  dry  density. 

Stress  dependence 

Figure  34  presents  a  set  of  logfm  vs  1  /a  curves  for  the  test  silt  with  medium  density  at  vari- 


ous  temperatures.  The  shapes  of  these  curves  are  similar  to  those  of  the  log£m  vs  1/a  curves 
shown  in  Figure  21,  and  hence  they  can  be  described  for  short-term  creep  by 


/m  =  t .  exp  A:, 


'(i  -  ±)| 


(20) 


where  t.  is  a  reference  constant  independent  of  temperature.  For  the  frozen  Fairbanks  silt 
with  densities  varying  from  1.18  to  1.23  g/cm\  t.  =  0.19  min.  The  parameter  k,  as  a  func¬ 
tion  of  temperature  is  shown  in  Figure  35  and  can  be  determined  by 


k,  =  49.7(0/0o)075 
k,  =  40.3(0/00)'  02 
For  long-term  creep, 

tm  =  tc  exp  k;  4-  - 


for  -0.5°  20>  -2°C 
for  -2°  2  0  >  -7°C. 


(22) 


(21) 


where  tc  is  the  time  to  failure  corresponding 
to  tc,  which  is  independent  of  temperature. 
From  the  test  data,  fc  has  a  magnitude  of 
about  900  min  for  Fairbanks  silt  with  medi¬ 
um  density.  The  parameter  k,'  depends  on 
temperature.  Its  variation  with  0  is  shown  in 
Figure  35  and  can  be  represented  by 


k!  =  1  18.8(0/0o)O9: 
for  -0.5°  2  0  >  -7°C. 


(23) 


10 


J _ i _ 


_i _ I _ 


-i _ I _ I  till 


0  I 


e 


/q  ,  Temperofure  Factor 


From  the  definition  of  creep  failure,  the 
creep  stress  a  in  eq  20  and  22  is  the  ultimate  Figure  35.  Parameters  k,  and  k,'  as  a  func- 
(maximum)  strength  (denoted  as  au|t)  of  fro-  tion  of  temperature. 
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Figure  36.  Comparison  of  the  computed  strength-lime  curves  with  test 
data. 


zen  soil  when  the  time  of  loading  is  the  same  as  the  time  of  creep  failure.  Then,  from  eq  20 
and  22  one  can  write  the  following  strength  relaxation  equations: 


ffult  = 


k,  a, 

a,  ln(t/t ,)  +  At, 


for  t  <  tc 


ffuit  - 


k,'oc 

Oq  ln(//fc)  +  k>' 


for  l  >  tc. 


These  two  equations  can  be  used  to  predict  the  ultimate  strength  of  the  frozen  silt  with  medi¬ 
um  density  at  any  given  time  of  loading.  The  curves  predicted  by  the  two  equations  agree 
well  with  the  test  data  (Fig.  36). 

It  should  be  especially  noted  that  one  cannot  predict  long-term  strength  by  simply  extrapo¬ 
lating  short-term  creep  data.  To  predict  long-term  strength,  of  course,  eq  25  should  be  used. 

Following  general  criteria,  if  we  assume  that  the  100-year  strength  is  the  limiting  long-term 

strength  <7|t,  then  it  can  be  calculated  by  eq  25  with  t  = 
100  years  =  5.25  x  10’  min.  The  predicted  values  of  aj, 
Table  7.  Predicted  values  o  o)t  as  a  function  of  temperature  for  frozen  Fairbanks  silt 
for  frozen  Fairbanks  silt  with  ..  .  . 

..  .  _  with  the  dry  densities  varying  from  1.18  to  1.23  g/cm* 

mediant  density  at  various  tem-  ,  ®.  . 

pcratures  are  shown  in  Table  7  and  Figure  37,  and  can  be  ex¬ 

pressed  as 


e 

CC) 

°u 

(kg/cm1) 

e 

CC) 

( kg/cm ') 

-0.3 

1.95 

-3 

8.81 

-1 

3.50 

-5 

14.22 

-2 

6.20 

-7 

18.67 

oh  =  3.49(0/0,)° 


where  o|t  is  in  kg/cm*. 

The  predicted  values  of  o\t  for  the  Fairbanks  silt 
were  compared  with  those  for  the  same  (or  similar) 


\*  *.7  ■.*  O  «.*  -  ’ 


Figure  37.  Predicted  limiting  long-term  strength  Figure  38.  Limiting  long-term  strengths  pre- 


as  a  function  of  temperature. 


dieted  by  various  investigators. 


material  with  different  water  contents  (Fig. 
38).  The  values  predicted  by  eq  25  are  con¬ 
siderably  higher  than  those  for  undisturbed 
ice-rich  frozen  silt  reported  by  McRoberts 
et  al.  (1978).  This  is  reasonable  because  the 
dry  densities  for  the  latter  are  much  higher 
than  those  for  the  former.  The  values  of  cti, 
predicted  by  eq  25  should  be  lower  than 
those  reported  by  Sayles  and  Haines  (1974), 
because  the  dry  densities  of  samples  pre¬ 
pared  in  his  investigation  are  higher  than 
those  in  this  study.  However,  Figure  38 
shows  that  the  results  are  almost  the  same. 
This  is  probably  because  the  predicted  val¬ 
ues  of  ff|t  in  Sayles’s  report  were  obtained 
by  extending  short-term  creep  data,  which 
may  considerably  underestimate  long-term 
strength. 


Temperature  dependence 
By  plotting  a  graph  of  log (tm/T)  against 
l/T  from  the  test  data  (Fig.  39),  it  was 

found  that  the  temperature  dependence  of  time  to  failure  tm  can  be  also  evaluated  in  terms  of 
an  exponential  expression: 


u 


tm  =  AfT  exp 


')] 

/?\7 

rJJ 

(27) 


where  A2  is  a  constant  equal  to  2.8  x  10'*  min/K  for  the  silt  with  medium  density. 
To  compute  Ue,  eq  27  is  rewritten  as 


,  Tim*  to  Foilur*  (min) 


where  Bt  =  -  log.4 ,  =  3.55.  The 
values  of  Ut  from  the  test  data  cal¬ 
culated  with  eq  28  are  plotted  in  Fig¬ 
ure  40  in  terms  of  log(/e  vs 
log(a/0,  + 1).  Obviously  the  data  in 
this  graph  can  be  also  fitted  by  eq  15 
very  well.  By  linear  regression  anal¬ 
ysis,  we  found  from  Figure  40  that 
AF  =  8100  kcal/mole  and  6  =  1.21, 
which  are  very  close  to  those  from 
Figure  28.  This  implies  that  there  is 
a  definite  relationship  between  mini¬ 
mum  creep  rate  and  time  to  failure. 
This  will  be  discussed  in  detail  in  the 
next  section. 


Figure  40.  Plot  of  log  U,  vs  logfo/o,  + 1). 


Effect  of  dry  density 
Figure  41  shows  time  to  failure  An 

as  a  function  of  dry  density  y d  for  different  stresses  at  -2°C.  Dry  density  clearly  has  a  signifi¬ 
cant  influence  on  w  The  denser  the  saturated  frozen  silt,  the  longer  the  time  to  failure.  Es¬ 
pecially  at  lower  stresses,  the  time  to  failure  for  the  dense  frozen  silt  is  much  longer  than  for 
ice-rich  frozen  silt. 

Figure  42  is  a  plot  of  logtm  vs  logfer/cr,)  for  the  low-density  samples  at  -2°C.  The  result  is  a 
straight  line,  so  the  time-to-railure  equation  for  low  density  has  the  form 


Figure  41.  Time  to  failure  as  a  function  of  dry  Figure  42.  Plot  of  logtm  vs  log  (ala,)  for  the 
density  for  different  stresses  at  -2°C.  specimens  with  low  density  at  -2°C. 
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where  is  the  time  to  failure  when  a  =  a, 
=  1  kg/cm1.  For  the  frozen  silt,  =  4.38  x 
10*  min  and  n  =  6.76  for  8  =  -  2°C. 

Figure  43  is  a  plot  of  log/m  vs  1/a  for  the 
high-density  samples  at  -2°C.  The  curve  in 
this  graph  seems  to  break  at  a  critical  stress 
ac  of  about  11.5  kg/cm2;  the  curve  can  be  de¬ 
scribed  by 

fm  =  /mo  tk,/°  (30) 

where  /mo  and  k„  depend  upon  the  stress 
level  for  a  given  temperature.  For  example, 
for  8  =  -  2°C,  tmo  =  0.03  min  and  k„  = 
116.4  kg/cm2  when  a  >  ac,  and  /mo  =  0.4 
min  and  Ar0  =  109.5  kg/cm2  when  a  <  ac. 


Figure  43.  Plot  of  logtm  vs  1/a  for  the 
specimens  with  high  density  at  -2  °C. 


Relationship  between  and  /„ 

In  each  plot  in  Figure  17,  one  can  draw  a  straight  line  that  nearly  passes  through  all  the 
minimum  points  of  logi  vs  log/  curves  over  a  certain  range  of  testing  stresses.  Similar  results 
were  reported  by  Martin  et  al.  (1981)  for  frozen  sand.  This  suggests  that  there  is  a  unique  re¬ 
lationship  between  the  minimum  creep  rate  and  the  time  to  failure  tm  over  a  certain  range 
of  stresses.  A  plot  of  all  minimum  points  for  various  temperatures  for  medium-density  sam¬ 
ples  (Fig.  44)  shows  that  this  relationship  is  identical  for  all  test  temperatures  and  can  be  ex¬ 
pressed  as 


im  =  0.086/m06 


(31) 


tm.  Time  to  Failure  (min) 


Figure  44.  Plot  of  logtm  vs  log tm  for  the  speci¬ 
mens  with  medium  density. 


where  tm  is  in  minutes,  and  im  is  in  min'1.  From  Figure  17i,  this  relationship  for  low-density 
samples  at  -2°C  is 

fm  =  0.0083 r£0i.  (32) 

In  Figure  17h,  one  cannot  draw  a  straight  line  through  all  of  the  minimum  points  over  the 
whole  range  of  test  stresses  because  of  the  large  drop  of  tm  at  lower  stresses.  However,  if  on¬ 
ly  a  certain  range  of  stresses  (say,  a  >  ac )  is  considered,  the  unique  relationship  between  £m 
and  /m  may  still  exist  and  can  be  written  as 


0.18^“. 


(33) 


Thus,  the  relationship  between  <?m  and  tm  has  the  general  form 


or 


=  C 


(34) 


where  p  depends  only  on  dry  density,  and  C  depends  on  the 
unit  of  time  and  the  dry  density  if  p  is  not  1 .  The  values  of  p 
and  C  for  different  ranges  of  dry  density  are  summarized  in 
Table  8. 

Table  8  shows  that  the  value  of  p  for  the  lower  densities  is 
very  close  to  1,  the  same  as  the  reported  value  for  polycrys¬ 
talline  ice  (Mellor  and  Cole  1982,  Cox  et  al.  1984).  However, 
it  is  significantly  greater  than  1  for  the  highest  density. 
Therefore,  the  simple  form  tm  <x  \/tm  is  not  applicable  to 
dense  frozen  soil. 


Table  8.  Values  of  p  and 
C  in  eq  34  for  different 
dry  densities.  The  unit  of 
time  is  minutes. 


id 


(g/cm') 

P 

C 

1.07-1.10 

1.01 

0.0083 

1.18-1.23 

1.06 

0.0860 

1.38-1.41 

1.16 

0.1800 

Creep  failure  strain  and  failure  criterion 

The  creep  failure  strains  «f  of  samples  with  medium  density  as  a  function  of  for  various 
temperatures  are  plotted  in  Figure  45.  Within  the  accuracy  of  the  test,  it  seems  that  the  fail¬ 
ure  strain  for  medium-density  samples  does  not  depend  on  and  0  over  a  wide  range  of  £m. 
However,  some  samples  did  not  fail  until  the  strain  was  more  than  15Vo  for  high  tempera¬ 
tures  ( a  -1°C)  and  low  stresses. 

Figure  46  shows  the  creep  failure  strains  as  a  function  of  for  different  dry  densities  at 
-2°C.  It  clearly  shows  that  the  failure  strain  strongly  depends  on  dry  density.  For  example. 


IO'8  IO'7  IO'6  to"5  I0'4  1C5  IO'2 


Em  ,  Minimum  Slroin  Role  (s') 


Figure  45.  Failure  strain  as  a  function  of  tmfor  specimens  with  medium  density  at 
various  temperatures. 
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Figure  46.  Failure  strain  as  a  function  of  tm  for  specimens  with 
various  dry  densities  at  -2  °C. 


Table  9.  Average  val¬ 
ues  of  6f  for  frozen 
Fairbanks  silt  with 
various  dry  densities 
at -2°C. 

( g/cm V  </ 


1.07-1.10 

1.18-1.23 

1.38-1.41 


0.0101 

0.0870 

0.1760 


•From  the  constant-stress 
creep  tests. 

the  failure  strain  decreases  from  about 
18%  to  1%  as  dry  density  decreases  from 
1 .4  to  1 .08  g/cmJ  at  the  higher  range  of 
#m.  Therefore,  it  is  incorrect  to  take  20% 
strain  as  a  creep  failure  criterion  for  plas¬ 
tic  frozen  soils  without  taking  into  consid¬ 
eration  the  density  or  water  content.  Dif¬ 
ferent  failure  criteria  for  frozen  soil  should 
be  selected  in  accordance  with  the  density 
and  the  failure  mode. 

Figure  46  shows  that  the  creep  failure  strain  for  low-density  samples  does  not  vary  with  im. 
It  has  an  average  value  of  about  1%,  which  is  very  close  to  that  of  polycrystalline  ice  (Mellor 
and  Cole  1982).  But  the  failure  strain  for  high  densities  varies  with  4m;  ef  suddenly  decreases 
from  about  16%  to  10%  as  tm  decreases  from  2.2  x  10'*  to  5.0  x  10'*  s'1.  When  <!m  is  greater 
than  2.2x10"*  s'1,  tf  increases  slightly  with  increasing  but  can  be  considered  practically 
constant.  Its  average  value,  tc^ther  with  that  of  samples  of  medium  and  low  density  at 
-1°C,  is  illustrated  in  Table  9  and  Figure  47.  From  Figure  47  the  creep  failure  strain  of  sat¬ 
urated  frozen  silt  as  a  function  of  dry  density  can  be  simply  evaluated  by 


Figure  47.  Failure  strain  as  a  function  of  dry 
density. 


ef  =  54.5(yd  -  7<i)  +  1 


(35) 
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Table  10.  Average 
values  of  tf-e»  for 
different  dry  densi¬ 
ty  groups  at  -2°C. 

yd 

(%/cm') _ irt. 


1.07-1.10  0.0086 

1.18-1.23  0.084 

1.38-1.41  0.173 


where  is  in  percent,  and  is  the  dry  density  at  which  tf  =  1%. 
For  the  frozen  silt  tested,  7<j  =  1.08  g/cm\ 

Using  the  values  of  the  initial  tangent  modulus  £;  of  the  frozen 
silt  for  constant-strain-rate  tests  shown  in  Table  4,  we  calculated 
the  instantaneous  strain  t0  by  simply  using  «0  =  o/E\.  Then  the 
quantity  ef-  e0  can  be  computed  for  each  specimen  (Table  10).  The 
values  of  ef-e0  are  nearly  equal  to  the  values  of  Cin  Table  8.  Thus, 
substituting  t(-t„  for  C  in  eq  34  we  obtained  the  following  creep 
failure  criterion: 


fm  ~  ef  fo- 


(36) 


Neglecting  e0,  eq  36  becomes 

*n.4=<f-  (37) 

Again,  the  unit  of  time  in  this  equation  must  be  in  minutes  if  p  is  not  equal  to  1 .  This  criter¬ 
ion  is  very  useful  because  it  links  three  basic  creep  failure  parameters  in  a  simple  mathemati¬ 
cal  form.  When  p  is  1,  eq  37  reduces  to 

*m  =  (38) 

which  is  exactly  the  same  as  Assur’s  (1980)  criterion  for  ice. 

Ladanyi  (1972)  proposed  a  criterion  for  frozen  soil  of  the  same  simple  form  as  eq  38.  As 
discussed  above,  applying  eq  38  to  ice-rich  frozen  soil  is  acceptable,  but  applying  it  to  dense 
frozen  soil  will  lead  to  significant  error  . 

Creep  model  and  prediction  of  creep  strain 

To  predict  the  time-dependent  deformation  of  frozen  soil,  a  number  of  researchers  have 
been  working  for  many  years  to  establish  a  quantitative  creep  model.  Since  the  creep  process 
of  frozen  soil  is  very  complicated,  especially  for  warm  plastic  frozen  soil,  it  is  difficult  to  de¬ 
velop  a  mechanistically  based  quantitative  model  involving  all  the  primary  creep  variables.  It 
is  possible,  but  the  model  must  be  exceedingly  complex.  Therefore,  empirical  modeling  is  by 
far  the  most  popular  approach. 

To  predict  creep  deformation  for  frozen  soil,  various  creep  models  have  been  developed. 
In  1962,  Vialov  proposed  a  primary  creep  model  that  has  been  widely  used  in  practice  and 
has  been  successfully  applied  to  various  frozen  soils  (Sayles  1968,  Sayles  and  Haines  1974). 
Ten  years  later,  Ladanyi  (1972)  developed  a  secondary  creep  model  (an  engineering  creep 
theory)  based  on  Hult’s  (1966)  creep  theory  on  metals.  Ladanyi’s  model  has  also  been  used  in 
engineering  because  of  its  simplicity. 

However,  great  care  must  be  taken  in  using  these  primary  and  secondary  creep  models, 
since  they  apply  only  for  the  duration  for  which  they  were  formulated.  For  example,  the  sec¬ 
ondary  creep  model  overestimates  creep  strain  during  primary  creep,  while  both  the  primary 
and  secondary  creep  models  underestimate  strain  during  tertiary  creep. 

A  similar  tertiary  creep  model,  developed  from  three  different  perspectives  by  Assur 
(1980),  Fish  (1980)  and  Ting  (1981),  can  be  used  to  describe  the  entire  process  of  creep.  In 
this  study,  Assur’s  model  is  used  to  fit  the  creep  data.  Based  on  the  study  of  creep  of  ice, 
Assur  (1980)  developed  a  simple  creep  model  that  is  derived  from  a  differential  equation 
quantifying  the  physical  basis  of  creep  and  has  the  form 

HO  =  (W^Om/O0  e0l/'m 
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(39) 


Table  11.  Values  of  0  in  eq 
39. 


*• 


(g/cm’J 


Values  of  0 


a  >  oc  a  s  ac 


1.38-1.41  0.7 

1.18-1.23  0.23/£ 

1,07-1.10  0.3  £ 

•  (tm  <  1000  min). 


2  3  4  5 

T- inT-l 


Figure  48.  Determination  of  parameter  0. 


where  0  can  be  determined  by  taking  logarithms  of  both  sides  of  eq  39  and  plotting  In l  vs 
(r-lnr-  1),  as  shown  in  Figure  48,  in  which  t  =  t/tm.  We  found  in  this  investigation  that  0 
varies  with  time  to  failure  rm  and  dry  density  74.  Its  average  values  are  shown  in  Table  1 1 . 

The  curves  predicted  by  eq  39  were  compared  with  the  observed  log*  vs  log*  curves  as 
shown  in  Figures  17e,  h  and  i.  The  comparison  shows  that  eq  39  can  be  used  to  fit  the  test 
curves  very  well  for  low-density  samples,  as  well  as  for  medium-  and  high-density  samples 
for  short-term  creep,  but  it  does  not  fit  the  test  curves  as  well  for  long-term  creep.  As  is 
shown  in  Figures  17e  and  h,  it  considerably  underestimates  t  during  primary  and  tertiary 
long-term  creep,  because  of  the  large  drop  of  *m. 

Integrating  eq  39  and  taking  into  consideration  eq  37  and  instantaneous  strain  6,,  we  de¬ 
rived  the  following  creep  equation: 


3t/tm  (01/tm)1 

■0X2-0 )  +  (l-0)(2-0X3-0) 


-  ...].  (40) 


Equation  40  is  not  applicable  for  long-term  creep  of  frozen  silt  with  a  high  density  because 
of  the  invalidity  of  eq  39.  In  this  case  the  creep  equation  obtained  from  integrating  eq  39  has 
the  form 


0X2-0) 


(0t/tmy 

(1-0X2-0X3-0) 


All  parameters  in  eq  40  and  41  can  be  determined  by  using  the  equations  discussed  earlier  in 
accordance  with  the  stress,  temperature  and  dry  density  of  concern. 

A  comparison  was  made  between  the  computed  curves  by  eq  40  and  41  and  test  data  in 
Figure  49,  showing  that  eq  40  can  adequately  describe  the  entire  creep  curves  for  the  samples 
with  low  density,  as  well  as  for  short-term  creep  of  the  samples  with  medium  and  high  densi¬ 
ty.  But  eq  40  and  41  do  not  predict  long-term  creep  as  well.  They  considerably  underestimate 
creep  strain  during  primary  and  tertiary  creep.  As  a  first  approximation,  they  are  still  applic¬ 
able  for  describing  the  long-term  creep  deformation  of  frozen  soil. 
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t,  Tim«  (min) 


a.  Short-term  creep  for  medium-density  samples  at 
-0.5°C. 


t,Tlm«  (min) 

b.  Short-term  creep  for  medium  density  samples  at 
-1  °C. 


c.  Short-term  creep  for  medium-density  samples  at  d.  Short-term  creep  for  medium-density  samples  at 
-2°C.  -3°C. 


e.  Short-term  creep  for  medium-density  samples  at  -5  °C. 
Figure  49.  Comparison  between  calculated  and  test  curves. 
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/.  Short-term  creep  for  medium-density  samples  at  -7°C. 


-10°C. 


Figure  49  (cont’d).  Comparison  between  calculated  and  test  curves. 


E,  Tru«  Strain 


k.  Long-term  creep  at  -2  °C  for:  a)  yd  =  1.230  g/cm1 ,  a  =  7.9  kg/ 
cm1;  b)  yd  =  1.185  g/cm1,  a  =  7.6  kg/cm1. 


I.  Long-term  creep  at  -7°C.  m.  Long-term  creep  for:  a)  yd  =  1.410  g/cm1,  a  =  11 

kg/cm1;  b)  yd  ~  1.407  g/cm1,  a  =  10  kg/cm1. 

Figure  49  (cont’d). 


STRENGTH  BEHAVIOR 

General  stress-strain  behavior  and  failure  mode 

The  constant-strain-rate  uniaxial  compression  tests  show  that  none  of  the  specimens  failed 
in  a  typically  brittle  mode,  but  they  failed  in  a  ductile  mode  at  the  testing  ranges  of  tempera¬ 
ture  and  strain  rate.  That  is,  fracturing  did  not  occur  until  the  axial  strain  approached  more 
than  25%.  However,  the  results  show  that  there  is  no  sharp  difference  between  ductile  and 
brittle  failures.  In  general,  the  failure  mode  gradually  shifts  from  one  to  the  other  as  the  test 
conditions  (such  as  temperature,  strain  rate  and  density)  change.  This  transition  can  be  clear¬ 
ly  seen  by  comparing  a  series  of  stress-strain  curves.  For  example,  from  Figure  14,  which 
presents  a  set  of  stress-strain  curves  for  specimens  with  medium  density  for  various  strain 
rates  at  -3°C,  one  can  see  that  the  brittleness  of  the  samples  gradually  increases  with  increas¬ 
ing  strain  rates.  A  more  significant  change  in  brittleness  seems  to  occur  at  a  strain  rate  above 
1.1  x  ]0~‘  s'1.  As  strain  rates  become  greater  than  1.1  x  10'5  s'1,  the  stress-strain  curves  drop 
sharply  after  reaching  a  peak,  indicating  a  great  increase  in  the  brittleness  of  the  soil. 
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Figure  50.  Typical  stress-strain  curves  Figure  51.  Typical  stress-strain  curves  for  various  dry 
for  various  temperatures  at  a  strain  rate  densities  at  a  strain  rate  of  1.1  x  10  '  s~'  and  a  tempera- 
of  1.1x10  ' s  '  for  medium-density  sam-  lure  of  -2  °C. 
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a.  Cracks  in  a  low-density  sample.  So  visible  cracks  in  a  medium-density 

sample. 

Figure  52.  Fairbanks  silt  specimens  with  different  densities  after  testing  at  a  fast  machine  speed  of  5 
cm/ min  and  a  temperature  of  -2  °C. 
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Figure  53.  Fairbanks  silt  specimens  of  medium  density  after  testing  at  a  slow  speed  of  0.001  cm/ min., 
showing  x  -shaped  shear  cracks. 

Figure  50  illustrates  a  set  of  stress-strain  curves  for  the  samples  with  medium  density  for 
various  temperatures  at  a  strain  rate  of  1.1  x  10'5  s'1.  It  shows  that  the  samples  behave  as  if 
the  soil  were  more  brittle  as  the  temperature  decreases.  The  tests  show  that  the  changes  in 
density  may  also  cause  a  transition  of  the  failure  mode.  Figure  5 1  shows  a  typical  set  of  stress- 
strain  curves  for  three  dry  densities  at  a  strain  rate  of  1.1  x  10'’  s'1  and  a  temperature  of 
-2°C.  It  is  obvious  that  the  lower  the  dry  density,  the  more  brittle  the  soil.  Figure  52  provides 
clear  evidence  to  confirm  this  conclusion.  Some  visible  cracks  can  be  seen  on  the  surface  of 
the  samples  with  low  density  (Fig.  52a),  but  no  cracks  can  be  seen  on  the  samples  with  high 
density  (Fig.  52b)  at  the  same  test  temperature  and  strain  rate. 

The  occurrence  of  cracks  in  a  specimen  does  not  always  indicate  that  it  failed  in  a  brittle 
mode.  In  fact,  as  discussed  earlier,  ductile  failure  sometimes  may  also  accompany  the  forma¬ 
tion  of  larger  visible  cracks.  For  instance,  the  large  x -shaped  cracks  seen  in  Figure  53 
resulted  from  a  typical  plastic  failure  of  the  samples  deformed  under  a  slow  strain  rate  of 
about  1.1  x  10'6  s  '. 

Peak  compressive  strength 

Figure  54  shows  a  plot  of  logam  vs  log(0/0o)  for  various  applied  strain  rates  e.  It  is  clear 
from  this  figure  that  the  peak  compressive  strength  am  of  the  frozen  silt  significantly  increas¬ 
es  with  decreasing  temperature  and  increasing  strain  rate.  As  reported  by  Sayles  and  Haines 
(1974),  the  peak  strength  of  frozen  soil  am  as  a  function  of  temperature  can  be  written  as 

om  =  A(»/Ho)m  (42) 

where  A  is  an  empirical  parameter  with  the  dimension  of  stress,  and  m  is  a  dimensionless 
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Figure  55.  Parameters  m  and  I/A  as  a  function  of  strain  rate. 


Tabic  12.  Values  of  A  and  m  in  eq 
43  for  frozen  Fairbanks  silt  with 
y4  *  1.18-1.23  g/cm1  at  6  a  -7°c. 


parameter.  For  a  particular  soil  type  and  a  certain  range 
of  temperature,  parameters  A  and  m  depend  on  strain 
rates.  The  values  of  A  and  m,  which  were  obtained  by 
linear  regression  analysis,  for  various  strain  rates  at  0 
2.  -7°C  are  listed  in  Table  12  (including  the  values  of 
correlation  coefficient  r)  and  plotted  in  Figure  55. 

There  is  a  significant  change  in  the  functional  rela¬ 
tionships  of  A  and  m  with  the  strain  rate  at  i  = 
1.1  xlO"5  s'*.  This  may  imply  that  different  deforma¬ 
tion  mechanisms  are  dominant  at  different  ranges  of 


strain  rate.  For  4  >  1.1  x  10"J  s'1,  the  deformation  process  of  the  frozen  silt  behaves  as  if  the 
soil  were  more  brittle,  while  for  i  <  1.1  x  lO-2  s'1,  it  behaves  as  if  the  soil  were  more  ductile. 
Hawkes  and  Mellor  (1972)  reported  similar  observations  for  fine-grained  polycrystalline  ice. 
Therefore,  we  suggest  that  a  strain  rate  of  about  1.1  x  10'1  s'*  be  considered  the  critical  strain 
rate  characterizing  the  ductile-brittle  transition  of  the  failure  mode  for  frozen  silt. 

From  Figure  55,  parameter  A  as  a  function  of  strain  rate  e  can  be  expressed  as 


where  h,  =  1.1  x  10"’  s"‘  is  the  strain  rate  hypothesized  to  be  the  transition  of  failure  mode 
from  ductile  to  moderate  brittle  for  this  frozen  silt,  A0  =  19. 1  kg/cm2  is  the  value  of  A  when 
l  =  and  a  is  the  reciprocal  of  the  slope  of  the  Ine  vs  l/A  curve  and  has  a  value  of  40 
kg/cm2  when  1.1  xlO"6  s'1  <  t  <  to  and  86  kg/cm2  for  i  >  e„. 

Parameter  m  as  a  function  of  k  can  be  written  as 

m  =  m0  +0.0517  ln(£0/£)  for  1.1  xl0"‘  s'1  <  e  <  e0 


m  =  m„  =  0.49 


for  t  >  l0 


Combining  eq  42  and  43  yields  a  constitutive  equation  as  follows: 

_  c‘Ao(e/0o)m 

m  a +/U  ln(£0/£) 

From  eq  4,  which  is  obtained  from  the  constant-stress  c  eep  test  data,  a  similar  constitutive 
equation  can  be  derived. 

It  is  obvious  from  Figure  54  that  the  peak  strength 

- 1-^— r — |  t  |  i  1 1 1 - n  increases  more  rapidly  as  the  temperature  decreases 

\  below  -7°C.  From  Figure  Al,  Appendix  A,  the  un- 

inri _  .  x.  .  _  frozen  water  content  W„  (%)  of  the  samples  with  me- 

_  '  Xt'62nl0  S  - 

—  '*  xJ  -  dium  dry  density  (corresponding  to  a  water  content 

I  -  \  V.  -  of  40.5%)  for  various  temperatures  can  be  calculated 

f  1  by 

I  XXXr  _  K  =  6.90(0/00)-° 47.  (46) 


y *1.18*1.23  g/cm3 

'd 


,  Unfrozen  Water  Content  (%) 

Figure  56.  Log  om  vs  logWu  for 
various  strain  rates. 


Equation  46  was  used  to  calculate  values  of  for 
the  various  testing  temperatures  shown  in  Figure  54 
to  produce  Figure  56,  using  the  same  peak  strength 
data.  All  of  the  curves  in  this  graph  also  break  at  a 
Wu  corresponding  to  -7°C.  If  the  peak  strength  is 
evaluated  in  terms  of  unfrozen  water  content,  a  func¬ 
tional  relationship  similar  to  that  relating  the  peak 
strength  to  temperature  can  be  easily  obtained. 

The  variation  of  peak  strength  with  dry  density  for 
various  strain  rates  at  -2°C  is  shown  in  Figure  57. 
Under  a  relatively  high  strain  rate,  the  peak  strength 
remains  almost  constant  over  the  test  range  of  dry 
density  of  1.07-1.43  g/cm’  at  -2°C.  But  under  rela- 
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Figure  57.  Peak  strength  as  a  function  of  dry  density  for  various 
strain  rates  at  -2  °C. 


tively  low  strain  rates,  the  peak  strength  signifi¬ 
cantly  increases  with  increasing  dry  density  for  the 
higher  dry  density  range.  For  lower  strain  rates, 
there  is  enough  time  for  the  sample  to  mobilize  its 
interparticle  frictional  resistance.  The  higher  the 
dry  density,  the  greater  the  frictional  resistance,  so 
the  higher  the  peak  strength. 

Table  13  compares  the  peak  strength  of  saturat¬ 
ed  and  partially  saturated  samples  of  medium  dry 
density.  The  peak  strength  of  the  partially  saturat¬ 
ed  samples  is  about  7-20%  less  than  that  of  satu¬ 
rated  ones,  even  though  the  saturation  degree  (in¬ 
cluding  ice  and  unfrozen  water)  for  the  partially 
saturated  soil  is  only  3-4%  less  than  for  the  satu¬ 
rated  soil.  This  indicates  that  for  a  given  density 
the  saturation  degree  may  significantly  influence 
the  peak  strength  of  frozen  silt. 


Table  13.  Peak  strength  of  saturated 
and  partially  saturated  samples  with 
an  average  dry  density  of  1.2  g/cm\ 


9 

rc) 

€ 

(S-) 

°m  (kg/cm') 

Partially 
saturated • 

Saturated t 

-i 

1.2x10' 

16.4 

19.3 

-i 

1.2x10— 

9.4 

10.1 

-2 

1.2x10-' 

20.6 

25.6 

-2 

1.2x10* 

13.8 

14.8 

-3 

1.2x10- 

17.1 

19.7 

-3 

1.2x10- 

9.7 

11.0 

*  Samples  saturated  under  a  normal  air  pressure, 
with  an  average  saturation  degree  of  about 
96-97 V.. 

t  Average  of  a  number  of  tested  samples. 


Initial  yield  strength 

Figure  58  shows  the  relationship  between  the  initial  yield  strength  ay  and  the  peak  strength 
am  for  various  testing  temperatures  at  dry  densities  ranging  from  1.18  to  1.23  g/cm'.  For  a 
particular  dry  density  there  is  a  unique  relationship  between  cry  and  om,  which  for  this  investi¬ 
gation  does  not  depend  on  temperature  or  strain  rate.  By  linear  regression  analysis  on  the  test 
data,  it  can  be  simply  expressed  by  a  linear  equation: 


For  between  1.18  and  1.23  g/cm\  k  =  0.79  and  c  =  -0.6  kg/cmJ.  Since  the  small  value 
of  c  may  result  from  the  test  error,  it  can  be  neglected.  Thus,  eq  47  can  be  written  as 


Oy  —  k  (Tin*  (48) 

The  initial  yield  strength  vs  peak  strength  for  various  dry  densities  is  plotted  in  Figure  59. 
It  indicates  that  the  coefficient  k  in  eq  48  depends  significantly  on  dry  density.  Through  a 
multiple  linear  regression  analysis,  it  was  found  from  the  test  data  that  the  initial  yield 
strength  ay  of  frozen  silt  can  be  determined  in  terms  of  peak  strength  am  and  dry  density 
by  the  following  equation  (with  r 1  =  0.9865  for  69  data  points): 


0,-0  79  am 
(r*0  9939) 


v,*  118*123  9  /cm5 


Om  ,  Peak  Strength  (kg/cm2) 

Figure  58.  Relationship  between  am  and  ay  for  the  spec i 
mens  with  medium  density  at  various  temperatures. 
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Figure  59.  Relationship  between  am  and  ayfor  various  dry 
densities  at  -2  °C. 
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Figure  60.  Initial  yield  strength  as  a  function  of  dry 
density  for  various  strain  rates  at  -2°C. 
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Figure  61.  Initial  yield  strength  as  a  function  of  vo/u 
metric  ice  content  for  various  strain  rates  at  -2°C. 


where  ><j  is  in  g/cm1  and  <rm  is  in  kg/cm2.  Neglecting  the  third  term  on  the  right  side  of  eq  49, 
we  have  the  simple  form  of  the  equation: 

<Ty  —  2.15am  -  1 . 1 5-yd  <rm  (50) 

or 

Oy  =  (2.15  —  1.157<i)<T|jj.  (51) 

This  equation  indicates  that  the  initial  yield  strength  ay  increases  with  decreasing  dry  densi¬ 
ty  for  the  saturated  frozen  silt.  This  can  be  clearly  seen  from  Figure  60,  which  presents  the 
yield  strength  as  a  function  of  dry  density  for  various  strain  rates  at  -2°C.  The  initial  yield 
strength  also  increases  with  increasing  volumetric  ice  content,  as  shown  in  Figure  61.  Figures 
60  and  61  also  show  that  the  initial  yield  strength  is  much  less  sensitive  to  dry  density  and  ice 
content  at  low  strain  rates  than  at  high  strain  rates  at  the  test  temperature.  Again,  this  simply 
means  that  the  initial  yield  strength  relaxes  much  faster  for  ice-rich  frozen  soil  than  for  dense 
frozen  soil.  Sayles  and  Carbee  (1981)  reported  a  similar  observation  for  the  same  soil  at  a 
temperature  of -1.67°C  and  strain  rate  of  5  x  10'1  s'1.  They  interpreted  this  stress  as  being  the 
stress  at  which  the  initial  fracture  of  pore  ice  in  frozen  soil  occurs.  This  conclusion  was 
reached  by  referring  to  the  stress-strain  behavior  of  polycrystalline  ice  where  the  ice  frac¬ 
tured  at  nearly  the  same  strain  (Hawkes  and  Mellor  1972). 

Equation  50  also  implies  that  the  difference  between  ay  and  am  decreases  with  decreasing 
dry  density.  Figure  51  illustrates  this  trend.  The  peak  strength  is  almost  twice  the  yield 
strength  for  7d  =  1.39  g/cm1,  but  there  is  no  significant  difference  between  them  for  >d  = 
1.10  g/cm1.  For  higher  dry  densities  the  fully  mobilized  frictional  resistance  between  soil  par¬ 
ticles  at  a  large  failure  strain  accounts  for  a  major  part  of  the  peak  strength  of  the  frozen 
soil.  For  lower  dry  densities,  the  strength  of  frozen  soil  is  substantially  governed  by  the  ice 
matrix,  which  fails  at  a  small  strain  of  about  1%  soon  after  its  plastic  yielding  (Sayles  and 
Carbee  1981). 

Failure  strain 

Test  results  of  failure  strain  tf  at  various  temperatures  and  strain  rates  for  medium-density 
samples  are  shown  in  Figure  62.  Within  a  certain  range  of  strain  rate  the  failure  strain  does 
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Figure  62.  Failure  strain  as  a  function  of  strain  rate  for  medium-density 
samples  at  various  temperatures. 
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Figure  63.  Variation  of  failure  strain  with  strain  rate  for  various  dry  densities 
at -2  °C. 


not  depend  on  temperature,  but  it  varies  slightly  with  strain  rate  and  has  a  maximum  value  at 
a  moderate  strain  rate  of  about  10’4  s’1.  However,  at  relatively  high  and  low  strain  rates,  the 
failure  strain  apparently  does  depend  on  temperature.  The  decrease  in  failure  strain  at  higher 
strain  rates  and  lower  temperatures  indicates  the  transition  of  failure  mode  from  ductile  to 
brittle. 

Figure  63  demonstrates  the  effect  of  dry  density  on  failure  strain  at  -2°C.  Failure  strain 
decreases  greatly  with  decreasing  dry  density.  For  each  density  group,  the  failure  strain  is  ap¬ 
proximately  constant  within  the  strain  rate  range  of  about  10‘*  to  6xl0’3  s’1.  The  average 
values  of  failure  strain  for  various  dry  density  groups  are  listed  in  Table  14  and  plotted  in 
Figure  47,  showing  a  good  agreement  with  the  values  from  the  creep  tests  discussed  earlier. 


Initial  yield  strain 

The  initial  yield  strain  discussed  here  refers  to  the  strain  at  the  yield  stress  ay.  Its  variation 
with  strain  rate  at  various  test  temperatures  for  medium-density  samples  is  shown  in  Figure 
64.  It  appears  that  for  a  given  dry  density,  initial  yield  strain  ey  significantly  depends  on 
strain  rate  but  is  almost  independent  of  temperature.  There  is  a  break  point  in  the  ty  vs  In* 
curve:  the  initial  yield  strain  increases  exponentially  with  increasing  strain  rate  when  t  >  10’’ 
s'1,  but  it  remains  almost  the  same  when  t  <  10'5  s'1.  The  same  observation  was  reported  by 
Mellor  and  Cole  (1982)  for  polycrystalline  ice.  This  evidence  strong- 
Table  14  Avera  e  ly  supports  Sayles  and  Carbee’s  (1981)  and  our  conclusion  that  the 

values  of  *€f  for  fro-  initial  yielding  of  frozen  silt  can  be  substantially  attributed  to  the 

zen  Fairbanks  silt  initial  fracture  of  the  ice  matrix.  Based  on  Figure  64,  initial  yield 

with  different  dry  strain  as  a  function  of  strain  rate  for  dry  densities  ranging  from 

densities  at  -2°C.  1.18  to  1 .23  g/cm’  can  be  expressed  as 


T  d 

(g/cm') 


I .08-1.10 
1.20-1.23 
I. 39- 1.43 


€y  =  0.27 +  0.045  In(;/10'5)  when  (  a  10'*  s'1 
and  (52) 

ey  =  0.27  when  i  <  10"’  s'1 


•From  constant-strtun- 
r«te  tests. 


where  ty  is  in  percent. 

Figure  65  shows  the  variation  of  the  initial  yield  strain  ey  with 
strain  rate  for  different  dry  densities  at  -26C.  There  seems  no  sig- 


tri 


Figure  64.  Initial  yield  strain  as  a  function  of  strain  rate  at 
various  temperature  for  medium-density  samples. 


Figure  65.  Initial  yield  strain  as  a  function  of  strain  rate  for  various 
dry  densities  at  -2°C. 


nificant  difference  of  «y  between  medium  and  high  densities.  However,  the  values  of  ey  for 
low-density  samples  appear  to  be  slightly  less  than  those  described  by  eq  52.  Its  magnitude 
ranges  from  about  0.2%  to  0.5%  as  strain  rate  increases  from  8.1  xlO"’  to  5.6  xlO'5  s'1, 
which  is  very  close  to  that  for  polycrystalline  ice  (Hawkes  and  Mellor  1972,  Mellor  and  Cole 
1982). 

Initial  tangent  modulus 

The  initial  tangent  modulus  £j  as  a  function  of  temperature  at  various  strain  rates  is  shown 
in  Figure  66.  Using  a  multivariable  regression  analysis,  we  obtained  the  following  equation 
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Figure  66.  Log E,  vs  log(6/0a)  for  medium-density 
samples  at  various  strain  rates. 


to  describe  the  initial  tangent  modulus  £j  as  a  function  of  temperature  and  strain  rate  for  74 
ranging  from  1.18  to  1.23  g/cm1: 

£i  =  6.1  x  lO'fe/tf,)0651  (t/t,)00*6  (53) 

where  £j  is  in  kg/cm1,  and  £,  is  a  reference  strain  rate  taken  as  1  s'1.  The  correlation  coeffi¬ 
cient  is  equal  to  0.8899,  and  the  standard  error  is  2.1  x  10J  kg/cm2. 

It  is  obvious  from  eq  53  that  the  initial  tangent  modulus  is  much  less  sensitive  to  strain  rate 
than  to  temperature.  Similar  results  were  reported  by  Hawkes  and  Mellor  (1972)  for  poly- 
crystalline  ice  and  by  Haynes  et  al.  (1975)  for  frozen  silt.  It  is  therefore  reasonable  to  con¬ 
sider  that  for  a  particular  material  the  initial  tangent  modulus  depends  only  on  temperature 
and  is  independent  of  strain  rate.  Thus,  eq  53  can  be  approximately  written  as  (with  r  = 
0.8690  and  s  =  2.1  x  10*  kg/cm2  for  52  data  points): 

£j  =  4.0  x  lO1(0/0,)°  (54) 

Figure  66  shows  that  the  data  on  E\  are  rather  scattered.  This  scattering  is  believed  to  be 
due  to  the  graphical  determination  of  E\  and  the  lack  of  stress-strain  data  at  the  initial  test 
period  for  higher  strain  rates.  It  is,  therefore,  essential  to  develop  standard  techniques  for 
determining  the  initial  tangent  modulus  and  for  precisely  recording  more  stress-strain  data 
for  the  fast  constant-strain-rate  tests. 

The  test  results  in  Table  4  show  that  there  is  no  significant  difference  in  E\  between  sam¬ 
ples  with  medium  and  high  densities  for  t  >  10"’  s'1  and  9  =  -2°C.  But  the  values  of  E\  for 
low-density  samples  are  considerably  greater  than  those  for  medium-density  samples.  E\  in¬ 
creases  from  about  5  x  101  to  2  x  104  kg/cm2  as  the  dry  density  decreases  from  about  1 .23  to 
1 .08  g/cm’  at  -2°C.  This  indicates  that  the  initial  tangent  modulus  of  dense  frozen  silt  can  be 
considerably  less  than  that  of  ice.  According  to  Hawkes  and  Mellor  (1972)  the  initial  tangent 
modulus  of  fine-grained  polycrystalline  ice  has  a  magnitude  of  2.8  x  104  to  7.4  x  104  kg/cm2 
at  -7°C.  We  have  shown  that  the  average  value  of  E\  for  frozen  Fairbanks  silt  with  an  aver- 
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age  dry  density  of  1 .20  g/cm3  ranges  from  1.3  x  104  to  2.2  x  104  kg/cm 2  at  the  same  tempera¬ 
ture. 

Plas-Tech  Corporation  (1966)  reported  that  the  modulus  values  of  Hanover  silt  at  -9.4°C 
are  from  6.2  x  10*  to  1.7  x  10*  kg/cm2,  which  is  close  to  that  of  the  tested  Fairbanks  silt.  Our 
test  results  of  the  initial  tangent  modulus  are  also  similar  to  those  reported  by  Haynes  et  al. 
(1973)  under  comparable  conditions. 

50%  peak  strength  modulus 

The  variation  of  50%  peak  strength  modulus  E,  with  temperature  for  various  strain  rates 
is  plotted  in  Figure  67  for  dry  densities  in  the  range  of  1 .20  to  1 .23  g/cm3.  The  50%  strength 
modulus  varies  with  both  temperature  and  strain  rate,  and  it  can  be  evaluated  by 

£,  =  ait/t^ie/do)”'.  (55) 

A  multiple  linear  regression  analysis  showed  that  the  parameters  in  this  equation  have  the 
values  of  o  =  7.6  xlO3  kg/cm3,  b  =  0.122,  and  n,  =  0.624  with  r  =  0.9447  and  5  =  1.1  x  103 
kg/cm3  for  54  data  points. 

The  50%  strength  modulus  as  a  function  of  strain  rate  is  shown  in  Figure  68  for  different 
dry  densities  at  -2°C.  Dry  density  obviously  has  a  significant  influence  on  the  50%  strength 
modulus.  The  curves  in  this  plot  are  nearly  parallel  to  each  other,  so  parameter  b  in  eq  55  is 
nearly  independent  of  dry  density  at  a  temperature  of  -2°C.  The  coefficient  a  is  the  only  par¬ 
ameter  that  determines  the  modulus  dependence  on  dry  density  if  the  parameter  n  in  this 
equation  is  also  independent  of  dry  density.  Assuming  that  b  does  not  depend  on  dry  density 
for  all  test  temperatures,  we  determined  from  the  test  data  that  a  -  9.3  x  103  kg/cm3  when  yj 
ranged  from  1 .08  to  1 . 10  g/cm3,  and  a  -  1 .2  x  103  kg/cm3  when  y^  ranged  from  1 .39  to  1 .43 


9/60,  Temperoture  Factor 

Figure  67.  Log  E,  vs  logfd/O* )  for  medium-density 
samples  at  various  strain  rates. 


Figure  68.  Log  E,  vs  logtfor  various  dry  densities  at  -2°C. 


In  eq  54  and  55  the  values  of  the  exponent  n,  of  (0/0o)  are  nearly  equal.  If  n,  is  considered 
to  be  invariant,  then  combining  these  two  equations  results  in  the  following  approximate  ex¬ 
pression: 


£j  =  0.54  (56) 

which  can  be  used  to  estimate  the  initial  tangent  modulus  of  frozen  silt  according  to  its  50% 
strength  modulus.  This  relation  is  helpful  for  determining  E\  because  is  easier  to  deter¬ 
mine. 

CORRESPONDENCE  BETWEEN  CONSTANT-STRESS  TESTS  AND 
CONSTANT-STRAIN-RATE  TESTS 

Mellor  and  Cole  (1982)  and  Sego  and  Morgens  tern*  recently  summarized  some  possible 
points  of  correspondence  between  constant-stress  tests  (creep  tests)  and  constant-strain-rate 
tests  (strength  tests)  on  polycrystalline  ice.  Ladanyi  (1981)  reported  that  this  correspondence 
is  also  valid  for  frozen  soils.  Our  work  shows  that  there  is  a  correspondence  between  the  re¬ 
sults  of  these  two  types  of  tests  on  frozen  silt.  The  results  of  the  failure  strain  and  am  (or  a)  vs 
t  (or  *„,)  curve  from  these  two  types  of  tests  are  nearly  the  same.  The  former  can  be  clearly 
seen  in  Figure  47,  in  which  the  failure  strains  from  both  the  creep  and  the  strength  tests  show 
a  good  correspondence  when  plotted  vs  dry  density.  The  latter  is  shown  in  Figures  69  and  70, 
in  which  am  vs  4  from  the  strength  tests  and  a  vs  im  from  the  creep  tests  are  plotted  together 
in  a  semi-log  scale,  showing  a  very  good  agreement  within  a  wide  overlap  range  of  strain 
rate. 

Because  of  the  correspondence,  one  can  use  the  constant-stress  and  constant-strain-rate 
test  data  equivalently  to  evaluate  the  time-dependent  stress-strain  behavior  of  frozen  soils.  In 
view  of  the  difficulties  in  gaining  creep  test  data  at  very  high  stresses  and  strength  test  data  at 
very  low  strain  rates,  it  might  be  appropriate  to  run  creep  tests  at  relatively  low  stresses  and 
strength  tests  at  relatively  high  strain  rates,  and  then  combine  the  results  from  the  two  types 
of  tests  to  establish  a  “complete”  a-t  curve.  Figure  71  represents  such  a  generalized  a-i 
curve  in  a  semi-log  scale  for  the  frozen  silt  tested.  For  discussion  purposes,  this  curve  is  di- 

*  Personal  communication,  University  of  Alberta,  1982. 


v  a*,. 


54 


Applied  Strom  Rote  (£)  or  Minimum  Creep  Rote  (Em)  (s  ) 

Figure  69.  Logtm  or  logt  vs  a  or  am  for  the  frozen  silt 
with  medium  density  at  -2  °C. 


=  118-1  23  g/cm 
0  *  -  iO°C 
o  Creep  Test 
•  Strength  Test 


I08  I0‘7  I0‘6  IO"5  I0"4  I0'3  lCfz  10 

Applied  Strom  Rote  (6)  or  Minimum  Creep  Rote  (Em)  (s  ') 

Figure  70.  Loglm  or  logl  vs  a  or  am  curve  for  the  frozen 
silt  with  medium  density  at  -10°C. 

vided  into  three  sections  by  two  characteristic  strain  rates,  tc  and  to,  corresponding  to  differ¬ 
ent  deformation  mechanisms.  As  discussed  earlier,  different  constitutive  equations  should  be 
used  to  describe  the  stress-strain  behavior  of  frozen  soil  in  accordance  with  different 
stress-temperature  Helds. 

We  could  not  get  the  upper  limit  of  the  a-t  curve,  that  is,  the  limiting  instantaneous 
strength  (denoted  as  oj„  in  Fig.  71),  in  this  investigation  because  of  the  head-speed  limit  of 
the  test  machine.  Many  other  researchers  may  also  have  this  kind  of  problem.  Therefore,  it  is 


■yt.f' 
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appropriate  to  define  an  identical  instan¬ 
taneous  strength  such  that  it  corre¬ 
sponds  to  a  certain  strain  rate.  In  the 
present  study,  instantaneous  strength 
(denoted  as  a0)  is  defined  as  the  peak 
strength  for  which  £  =  10" 1  s'1.  Then, 
from  eq  26  and  42  we  obtain 

<r„  =  0.18(0/9')°™  a,  (57) 

which  can  be  used  for  predicting  the  lim¬ 
iting  long-term  strength  for  saturated 
frozen  Fairbanks  silt  with  dry  densities 
ranging  from  1.18  to  1.23  g/cmJ. 


Figure  71.  Generalized  curve  of  logt  vs  a. 


CONCLUSIONS 

The  constant-stress  and  constant-strain-rate  compression  tests  lead  to  the  following  con¬ 
clusions. 

There  exists  a  critical  creep  strength  ac,  which  can  be  used  to  classify  the  creep  of  frozen 
silt  into  two  basic  types:  short-term  creep  for  which  a  >  ac,  and  long-term  creep  for  which  a 
<  uc.  Short-term  creep  is  controlled  by  glide  creep,  while  long-term  creep  is  governed  by  dis¬ 
location  creep.  They  obey  different  creep  laws  and  thus  different  constitutive  equations. 

The  rate  process  theory  is  applicable  for  frozen  soil.  However,  the  linear  function,  Ue  = 
AF-po,  is  not  valid  for  the  lower  stress  range.  Instead,  a  power-law  function,  Ue  = 
A.F(1  +  a/<j|)'4,  works  very  well  over  a  wide  range  of  stresses.  A  free  activation  energy  AFof 
about  8600  kcal/mole  was  obtained  for  frozen  Fairbanks  silt  with  medium  dry  density, 
which  is  much  higher  than  that  reported  by  other  investigators. 

Dry  density  has  a  significant  effect  on  the  creep  behavior.  The  flow  law  of  ice-rich  frozen 
silt  is  similar  to  that  of  polycrystalline  ice  and  can  be  described  by  a  simple  power  law  over  a 
wide  range  of  creep  stresses.  For  the  dense  frozen  silt,  though,  the  closely  packed  soil  parti¬ 
cles  greatly  impede  the  creep  process  within  secondary  creep  at  lower  stresses,  so  its  creep  law 
cannot  be  described  by  a  single  power  law  over  a  wide  range  of  creep  stress  but  requires  more 
than  one  constitutive  equation. 

Different  types  of  creep  have  different  strength-loss  equations.  One  cannot  predict  long¬ 
term  strength  by  simply  extrapolating  the  short-term  creep  data.  The  100-year  strength  q,  of 
frozen  Fairoanks  silt  with  medium  density  can  be  determined  by 

q,  .  3.49(0/0„)087  (58) 


or 


<71,  =  0.18(0/0,)° 38  a,  (59) 

where  q,  is  in  kg/cm*  and  q>  is  the  peak  strength  of  the  silt  at  1  =  10"’  s  '. 

Failure  strain  «f  does  not  seem  to  be  sensitive  to  temperature  within  the  accuracy  of  the  test 
data,  and  it  also  is  not  sensitive  to  4m  for  a  particular  failure  mode  or  deformation  mechan¬ 
ism.  However,  it  is  very  sensitive  to  dry  density.  It  increases  proportionally  from  about  1% 
to  18^0  as  dry  density  increases  from  1.08  to  1.40  g/cmJ. 
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The  creep  failure  criterion  for  this  frozen  soil  can  take  the  general  form  of  tm  <%  =  ([, 
where  p  depends  only  on  dry  density,  and  tm  is  in  minutes  if  p  is  not  1 .  For  frozen  Fairbanks 
silt,  the  values  of  p  are  given  in  Table  7. 

Assur’s  (1980)  creep  model  for  ice  can  well  describe  the  complete  creep  curve  for  ice-rich 
frozen  silt  as  well  as  the  short-term  creep  curve  for  dense  frozen  silt,  but  it  cannot  describe 
long-term  creep  as  well.  It  considerably  underestimates  creep  deformation  during  primary 
and  tertiary  creep  for  long-term  creep. 

The  frozen  Fairbanks  silt  samples  generally  failed  in  a  ductile  mode  under  the  constant- 
strain-rate  test  conditions.  But  at  very  high  strain  rates  and  low  temperatures,  a  moderate 
brittle  failure  may  occur. 

The  compressive  peak  strength  am  of  the  frozen  silt  is  very  sensitive  to  temperature  and 
strain  rate.  Both  A  and  m  in  the  equation  am  =  A($/6„)m  depend  on  the  strain  rate  and  tem¬ 
perature  range.  It  has  been  found  that  for  0  a  -7°C  both  the  1/A  vs  Ine  and  the  m  vs  lne 
curves  break  at  a  strain  rate  of  1.1  xlO*  s'1,  perhaps  indicating  the  transition  of  failure  mode 
from  ductile  to  brittle. 

For  the  material  tested  at  -2°C,  peak  strength  was  found  to  be  not  sensitive  to  dry  density 
at  relatively  high  strain  rates,  but  it  decreases  considerably  with  decreasing  dry  density  at 
lower  strain  rates. 

The  initial  yield  strength  <ry  is  sensitive  not  only  to  temperature  and  strain  rate,  but  also  to 
dry  density,  indicating  that  it  is  essentially  related  to  the  initial  fracture  of  pore  ice  in  frozen 
soil.  There  is  a  definite  relationship  between  am  and  ay  for  the  silt  that  is  not  dependent  on 
temperature  and  strain  rate: 

<jy  =  (2.15-1.157d)am.  (60) 

For  a  certain  range  of  applied  strain  rates,  the  failure  strain  of  the  silt  does  not  vary  with 
temperature,  but  it  depends  strongly  on  dry  density.  For  the  saturated  frozen  Fairbanks  silt, 
failure  strain  decreases  from  18.7%  to  1%  as  dry  density  decreases  from  about  1.40  to  1.08 
g/cm*  over  a  range  of  applied  strain  rates  from  1.1  x  10'*  to  1.1  x  10"’  s'1. 

The  initial  yield  strain  of  the  frozen  silt  does  not  depend  on  temperature  but  increases  with 
increasing  strain  rate  when  t  a  10*5  s'1.  When  i  <  10'’  s~\  it  remains  almost  constant.  It 
decreases  slightly  with  decreasing  dry  density.  For  the  samples  with  low  dry  unit  weight,  it 
has  a  magnitude  of  0.2-0. 5%  over  a  strain  rate  range  of  1.1  x  10'>  to  1.1  x  10'*  s'1. 

The  initial  tangent  modulus  E\  of  the  frozen  silt  is  not  sensitive  to  strain  rate  in  the  range  of 
ratios  used,  but  it  varies  with  temperature.  The  value  of  £j  for  the  tested  silt  with  ranging 
from  1.18  to  1.23  g/cm’  can  be  determined  (in  kg/cm2)  by 

E{  =  4.0  x  lO*(0/0«)0636.  (61) 

The  50%  peak  strength  modulus  E,  of  the  silt  as  a  function  of  temperature,  strain  rate  and 
dry  density  can  be  expressed  as 

E,  =  a{t/tl)b(e/6,r  (62) 

where  b  =  0.122,  n  =  0.624,  and  a  depends  on  dry  density  as  follows: 

td  a 

(g/cm1)  (kg/cm2) 

1.08-1.10  9.3  xlO1 

1.20-1.23  7.6x10* 

1.39-1.43  1.2x10* 
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The  initial  tangent  modulus  E\  and  the  50%  strength  modulus  £,  can  be  related  to  each 
other  by 

£i  =  0.54a/*,r0122£i.  (63) 

There  is  a  correspondence  between  the  constant-stress  and  constant-strain-rate  tests  for 
the  frozen  silt  under  the  test  conditions  in  this  study.  In  particular,  the  average  failure  strains 
over  a  certain  range  of  t  (or  £m)  from  the  two  types  of  tests  are  approximately  the  same,  and 
the  crm-£  curves  from  the  strength  tests  and  the  o-£m  curves  from  the  creep  tests  are  in  good 
agreement. 

Creep  tests  at  a  relatively  low  stress  and  strength  tests  at  a  relatively  high  strain  rate  can  be 
used  to  yield  a  complete  a-i  curve  for  a  soil. 
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APPENDIX  A:  UNFROZEN  WATER  CONTENT  DATA 


Table  Al.  Unfrozen  water  content  data  for  the  samples  with  three  typical  water  contents.  The 
data  are  plotted  in  Figure  Al. 


w  - 

30. 3X* 

W  - 

40.5Xt 

W  - 

49.9*7 

Unfrozen 

Unfrozen 

Unfrozen 

Temperature 

water  content 

Temperature 

water  content 

Temperature 

water  conten 

e,*c 

W  ,  X 

0,*C 

W  ,  X 

9,*C 

W  ,  X 

u 

u 

u 

-15.55 

3.06 

-18.45 

1.90 

-18.50 

2.17 

-14.12 

3.18 

-15.44 

1.98 

-15.57 

2.21 

-10.85 

3.13 

-13.15 

2.06 

-13.39 

2.23 

-8.95 

3.48 

-11.06 

2.13 

-11.09 

2.26 

-6.86 

3.55 

-8.40 

2.53 

-8.48 

2.71 

-4.96 

4.13 

-7.36 

2.27 

-7.44 

2.80 

-3.86 

4.43 

-6.15 

2.80 

-6.09 

3.03 

-2.62 

4.90 

-4.43 

3.19 

4.51 

3.56 

-1.76 

5.55 

-3.62 

3.65 

-3.52 

3.86 

-1.48 

6.42 

-2.36 

4.31 

-2.44 

4.54 

-1.06 

7.83 

-1.40 

6.25 

-1.61 

5.57 

-0.88 

7.95 

-0.82 

6.78 

-0.82 

6.90 

-0.54 

8.12 

-0.30 

13.77 

-0.28 

13.85 

-0.51  8.53 
-0.38  9.64 
-0.33  11.86 
-0.28  16.31 


*  After  Tice  (personal  communication), 
t  After  Xu  (peraonal  conaunlcatlon). 


Figure  Al.  Unfrozen  water  content  as  a  function  of  temperature:  I )  w 
=  30.3%  (Tice,  pers.  comm.);  2)  w  =  40.5%,  and 3)  w  =  49.9%  (Xu, 
pen.  comm.). 
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APPENDIX  B:  PHYSICAL  PROPERTIES  OF  SAMPLES  TESTED 


Sample 

No. 

Water 
content 
v  (Z) 

Dry 

density 

Yd  (g/cm3) 

Volume  Ice 
contentt 

Saturation 
degree** 
a  (Z) 

Void 

ratio 

e 

Porosity 

n 

5-53 

44.2 

1.162 

0.508 

98.1 

1.307 

0.567 

5-54 

40.9 

1.200 

0.470 

96.0 

1.233 

0.552 

5-55 

43.0 

1.191 

0.494 

96.2 

1.251 

0.556 

5-56 

40.7 

1.205 

0.469 

96.2 

1.224 

0.550 

5-57 

42.0 

1.186 

0.454 

96.0 

1.261 

0.558 

5-58 

40.2 

1.212 

0.440 

95.7 

1.211 

0.548 

3-59 

41.9 

1.202 

0.459 

99.5 

1.230 

0.552 

4-60 

43.5 

1.182 

0.472 

99.6 

1.265 

0.558 

4-61 

42.4 

1.200 

0.465 

99.2 

1.234 

0.552 

5-63 

41.3 

1.192 

0.447 

95.5 

1.247 

0.555 

5-64 

41.6 

1.191 

0.451 

95.9 

1.252 

0.556 

5-65 

41.2 

1.196 

0.447 

95.5 

1.242 

0.554 

5-66 

40.7 

1.200 

0.442 

95.3 

1.232 

0.552 

5-67 

41.2 

1.196 

0.472 

96.1 

1.243 

0.554 

5-68 

41.2 

1.197 

0.473 

96.2 

1.238 

0.553 

4-69 

43.1 

1.182 

0.467 

98.5 

1.265 

0.559 

4-70 

43.8 

1.181 

0.500 

100.0 

1.269 

0.559 

4-71 

42.2 

1.203 

0.500 

99.8 

1.228 

0.551 

4-72 

43.4 

1.184 

0.507 

99.7 

1.262 

0.558 

4-73 

41.5 

1.214 

0.495 

99.7 

1.207 

0.547 

4-74 

42.0 

1.205 

0.498 

99.7 

1.225 

0.551 

4-75 

43.2 

1.189 

0.495 

99.8 

1.254 

0.556 

4-76 

41.6 

1.210 

0.495 

100.0 

1.215 

0.549 

4-77 

42.0 

1.195 

0.494 

98.0 

1.242 

0.554 

5-78 

41.2 

1.199 

0.473 

96.3 

1.237 

0.553 

5-79 

40.5 

1.205 

0.466 

95.7 

1.225 

0.550 

4-80 

41.8 

1.206 

0.496 

99.3 

1.220 

0.550 

4-81 

45.2 

1.173 

0.525 

100.0 

1.285 

0.562 

5-82 

41.9 

1.186 

0.418 

95.3 

1.261 

0.558 

4-83 

46.0 

1.144 

0.522 

99.5 

1.343 

0.573 

5-84 

40.9 

1.200 

0.410 

95.2 

1.233 

0.552 

5-85 

42.7 

1.175 

0.425 

95.7 

1.281 

0.562 

6-86 

40.0 

1.240 

0.412 

99.0 

1.162 

0.537 

6-87 

40.9 

1.227 

0.419 

99.8 

1.183 

0.542 

6-88 

41.6 

1.214 

0.424 

99.6 

1.206 

0.547 

6-89 

42.0 

1.210 

0.428 

99.5 

1.214 

0.548 

6-90 

41.1 

1.222 

0.456 

99.4 

1.193 

0.544 

6-91 

40.4 

1.229 

0.449 

99.0 

1.179 

0.541 

6-92 

40.8 

1.224 

0.452 

99.2 

1.188 

0.543 

6-93 

41.9 

1.206 

0.425 

98.6 

1.220 

0.549 

6-94 

41.3 

1.219 

0.457 

99.4 

1.198 

0.545 

6-95 

41.5 

1.214 

0.423 

99.0 

1.206 

0.547 

6-96 

40.2 

1.237 

0.475 

99.8 

1.166 

0.538 

6-97 

40.6 

1.230 

0.478 

99.7 

1.179 

0.541 

6-98 

40.1 

1.238 

0.412 

99.0 

1.164 

0.538 

6-99 

42.7 

1.203 

0.506 

100.0 

1.227 

0.551 

6-100 

41.8 

1.214 

0.499 

100. 0 

1.207 

0.547 

6-101 

40.9 

1.226 

0.504 

100.0 

1.186 

0.542 

7-102 

41.4 

1.213 

0.505 

99.6 

1.209 

0.547 

7-103 

42. 1 

1.205 

0.511 

99.8 

1.224 

0.550 

7-104 

41.5 

1.214 

0.513 

100.0 

1.207 

0.547 

7-105 

40.2 

1.222 

0.493 

98.4 

1.191 

0.544 

7-106 

41.0 

1.218 

0.508 

99.3 

1.201 

0.546 

7-107 

42.2 

1.197 

0.515 

99.2 

1.239 

0.553 

7-108 

40.6 

1.214 

0.495 

98.6 

1.207 

0.547 

7-109 

41.4 

1.190 

0.501 

99.1 

1.250 

0.556 

7-110 

40.9 

1.221 

0.501 

99.5 

1.194 

0.544 
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# 


Sample 

No. 

Water 
content 
w  (Z) 

Dry 

density 

Yd  (g/cm3) 

Volume  ice 
contentt 

Saturation 

degree** 

8  (Z) 

Void 

ratio 

e 

Poroeity 

n 

7-111 

40.5 

1.214 

0.500 

98.7 

1.205 

0.546 

7-112 

41.2 

1.216 

0.510 

99.4 

1.204 

0.546 

7-113 

40.6 

1.221 

0.509 

98.8 

1.195 

0.544 

7-114 

41.2 

1.208 

0.512 

98.3 

1.219 

0.549 

7-115 

40.2 

1.227 

0.507 

98.8 

1.183 

0.542 

7-116 

41.6 

1.182 

0.506 

98.8 

1.265 

0.559 

7-117 

41.0 

1.218 

0.513 

99.5 

1.201 

0.546 

7-118 

41.7 

1.210 

0.519 

99.9 

1.215 

0.549 

8-119 

45.9 

1.168 

0.555 

100.0 

1.295 

0.564 

8-120 

41.8 

1.210 

0.521 

100.0 

1.214 

0.548 

8-121 

42.6 

1.198 

0.526 

100.0 

1.236 

0.553 

8-122 

43.7 

1.179 

0.532 

100.0 

1.272 

0.560 

8-123 

42.3 

1.202 

0.524 

100.0 

1.229 

0.552 

8-124 

42.4 

1.202 

0.520 

100.0 

1.231 

0.552 

8-125 

42.3 

1.203 

0.519 

100.0 

1.227 

0.551 

8-126 

43.3 

1.174 

0.519 

99.8 

1.282 

0.562 

8-127 

41.4 

1.211 

0.510 

99.6 

1.211 

0.548 

8-128 

43.1 

1.187 

0.522 

99.6 

1.257 

0.557 

8-129 

40.7 

1.222 

0.506 

99.4 

1.190 

0.543 

8-130 

43.) 

1.184 

0.515 

99.4 

1.262 

0.558 

8-131 

42.5 

1.195 

0.512 

99.5 

1.241 

0.554 

8-132 

40.6 

1.229 

0.501 

100.0 

1.181 

0.541 

8-133 

42.0 

1.203 

0.508 

99.5 

1.226 

0.551 

8-134 

41.4 

1.221 

0.508 

99.7 

1.195 

0.544 

8-135 

42.0 

1.211 

0.512 

100.0 

1.212 

0.548 

9-136 

43.0 

1.190 

0.516 

99.8 

1.252 

0.556 

9-137 

43.5 

1.182 

0.531 

100.0 

1.267 

0.559 

9-138 

44.2 

1.168 

0.522 

99.2 

1.293 

0.564 

9-139 

44.6 

1.165 

0.532 

99.7 

1.301 

0.565 

9-140 

41.7 

1.206 

0.494 

100.0 

1.221 

0.550 

9-141 

44.2 

1.171 

0.501 

99.4 

1.289 

0.563 

9-142 

42.9 

1.187 

0.502 

100.0 

1.257 

0.557 

9-143 

42.0 

1.190 

0.492 

97.8 

1.250 

0.556 

9-144 

42.8 

1.195 

0.504 

100.0 

1.241 

0.554 

9-145 

44.2 

1.200 

0.513 

100.0 

1.039 

0.510 

9-146 

41.6 

1.200 

0.479 

97.3 

1.233 

0.552 

9-147 

44.0 

1.173 

0.499 

98.9 

1.283 

0.562 

9-148 

43.3 

1.182 

0.435 

99.2 

1.265 

0.558 

9-149 

43.0 

1.177 

0.429 

97.0 

1.275 

0.560 

9-150 

42.9 

1.205 

0.473 

100.0 

1.223 

0.550 

9-151 

43.6 

1.184 

0.440 

99.1 

1.263 

0.558 

9-152 

43.6 

1.179 

0.438 

98.3 

1.272 

0.560 

10-153 

45.2 

1.155 

0.449 

98.3 

1.321 

0.569 

10-154 

47.8 

1.123 

0.468 

99.0 

1.385 

0.581 

10-155 

44.4 

1.163 

0.522 

100.0 

1.304 

0.506 

10-156 

45.2 

1.155 

0.529 

99.6 

1.320 

0.569 

10-157 

49.6 

1.096 

0.533 

99.4 

1.444 

0.591 

10-158 

49.8 

1.091 

0.533 

99.0 

1.457 

0.593 

10-159 

49.8 

1.085 

0.530 

98.9 

1.472 

0.595 

10-160 

50.1 

1.088 

0.535 

99.1 

1.462 

0.594 

10-161 

48.8 

1.105 

0.528 

99.1 

1.424 

0.588 

10-162 

49.3 

1.099 

0.531 

99.1 

1.439 

0.590 

10-163 

50.8 

1.070 

0.534 

98.1 

1.505 

0.601 

10-164 

49.8 

1.086 

0.531 

98.2 

1.468 

0.595 

10-165 

** 

00 

1.104 

0.524 

98.4 

1.425 

0.588 

10-166 

51.2 

1.080 

0.544 

100.0 

1.480 

0.597 

65 


V' 


:»v  •»  .v 


t  The  volume  ice  content  iv  is  defined  as 


1 


v 


_ Volume  of  ice _ 

Total  volume  of  soil  mass 


and  calculated  by 


(w  - 


w  ) 
u 


G.  Y 

i  w 


where 

G^  »  specific  gravity  of  ice,  equal  to  0.917 
Yw  «  unit  weight  of  water,  taken  as  1  g/cm3 
w  -  water  content  (?) 

wu  ■  unfrozen  water  content  of  samples,  which  can  be  found  in 
Appendix  A. 

**  The  saturation  degree  8^  is  defined  as 

Volume  of  ice  and  unfrozen  water 
i  Volume  of  pore 

and  calculated  by 

100[w-(l-G1)wu]G 
81  G^  ~ 

where  G  is  the  specific  gravity  of  the  soil  particles,  equal  to  2.680,  and 
e  is  the  void  ratio  of  samples. 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 
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